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1. Introduction 

Certain quantum field theories are known to possess a continuous set of inequivalent ground 
states, or in other words, a quantum moduli space. Every point in the moduli space cor- 
responds to a vacuum state with zero energy. Turning on a non-zero temperature T will 
generically "lift" moduli space, leaving a much smaller set of thermal equilibrium states. For 
some range of temperatures, there may be a unique equilibrium state. For other tempera- 
tures, there may be multiple degenerate equilibrium states related by spontaneously broken 
global symmetries. 

One may define a thermal effective potential, or free energy functional, using the same 
coordinates which parametrize the zero temperature moduli space. This will turn the flat 
T = zero energy surface into a non-trivial T > free energy surface. Equilibrium states 
correspond to the global minima of this free energy surface. Computing this free energy 
surface in an interacting theory is, of course, non-trivial. Supersymmetry provides little 
help, since cancellations between bosonic and fermionic particles in virtual processes are 
spoiled by their different statistics at non-zero temperature. Nevertheless, two extreme limits 
are interesting and amenable to analytical calculation: arbitrarily low temperatures and 
asymptotically high temperatures. In the former case, one might expect the free energy 
surface to be a slight deformation away from the flat surface of moduli space. What does 
this lift look like, and where do the minima of the free energy surface lie? In the high 
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temperature regime, thermal fluctuations should have a disordering effect on the system. Are 
spontaneously broken global symmetries restored? At what temperature? 

In this work, we examine the effects of thermal fluctuations on the equilibrium properties 
and realization of global symmetries in the simplest asymptotically free supersymmetric gauge 
theory with a continuous moduli space, SU(2) N = 2 supersymmetric Yang-Mills theory. 
Much about this theory is known from the celebrated work of Seiberg and Witten [1, 2]. The 
following features make it an attractive model for our purposes: 

(i) The quantum theory has a continuous moduli space of vacua — it is a one-complex 
dimensional Kahler manifold parametrized by a single complex number, u. 

(ii) Each vacuum describes a Coulomb phase where the long distance dynamics is Abelian; 
there is no vacuum state with long distance non- Abelian dynamics (i.e., confinement). 

(iii) Generic ground states spontaneously break a discrete ^-symmetry. 

(iv) Asymptotic freedom guarantees that vacua in the neighborhood of infinity on moduli 
space have weakly-coupled descriptions in terms of the light elementary fields. 

(v) Two distinguished "singular" points in moduli space exist where extra massless states 
with spin < | appear. The corresponding particles are magnetically charged under 
the long-distance Abelian gauge group and may be interpreted as magnetic monopoles 
or dyons. For vacua in neighborhoods of these special points, a low energy effective 
description is strongly-coupled in terms of the elementary fields. However, a version of 
electric-magnetic duality provides a weakly- coupled formulation in terms of dual fields. 

The combination of asymptotic freedom and electric-magnetic duality enables one to use 
weak coupling methods to explore the dynamics both near and far from the singular points 
in moduli space. There is a dynamically generated mass scale A in the theory (analogous 
to Aqcd)- For temperatures much greater than A, the free energy may be computed as an 
asymptotic expansion in the small effective gauge coupling g 1 (T) . For temperatures much less 
than A, one may use appropriate low energy effective descriptions near infinity, or near the 
special points on moduli space, to compute the free energy as an expansion in the appropriate 
effective gauge coupling (either g 2 (u) or its magnetic dual). 

Thermal effects in SU(2) N = 2 gauge theory at low temperature have previously been 
studied by Wirstam [3]. In this work, it was asserted that the free energy density was 
locally minimized asymptotically far out on moduli space, and on circles of non-zero radius 
surrounding the singular points. It was not made clear which local minima represented the 
global minimum. The free energy surface found in Ref. [3] is depicted on the left side of 
Figure 1. In this figure, arrows depict directions of free energy decrease (i.e., minus the 
gradient). The picture implies non-monotonic behavior as one moves from a singular point 
to infinity, with a free energy barrier separating the large u domain from the region near 
the singular points, and some sort of instability at the singular points. Such features are 
unexpected and surprising. One puzzle is why the free energy surface slopes downward to 
infinity. Massive states get heavier as one moves further out on moduli space so, in the 
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Figure 1: Qualitative form of the free energy surface in SU(2) M =2 Yang-Mills theory at T -c A. 
Arrows indicate directions on moduli space for which free energy decreases. The singular points at 
±ito are represented by heavy dots. The large dotted circle at infinity serves to guide the eye. Left: 
Asserted behavior from Ref. [3] . The dashed circles surrounding the singular points at ±uo represent 
valleys of stable local minima. Right: Results of our analysis. The singular points are stable minima. 



absence of interactions, one would expect their contribution to the pressure to decrease (since 
the associated particle density falls exponentially due to Boltzmann suppression). The free 
energy density is minus the pressure, so the decoupling of massive states as one approaches 
the boundary of moduli space should lead to a rising free energy. Do interactions, in an 
asymptotically free theory, really change this simple behavior? 

A second puzzle concerns the circle of minima around each singular point. What physical 
mechanism leads to this? There is no continuous global symmetry whose action on moduli 
space produces phase rotations around a singular point, and whose spontaneous breaking 
could explain such a circle of free-energy minima. 

The main purpose of this paper is to derive the correct behavior of the free energy surface 
at low temperatures by reconsidering the computations of Ref. [3]. Using effective field theory 
techniques, we systematically evaluate the contributions of successively longer wavelength 
fluctuations to the effective scalar potential. The free energy surface, viewed as a functional 
of the translationally invariant expectation value which parameterizes moduli space, may be 
identified with this effective potential. Unlike Ref. [3], we find the simple behavior sketched 
on the right side of Figure 1. In the low temperature regime, the asymptotic region of the free 
energy surface is locally unstable; the free energy decreases as one moves inward from infinity. 
The two singular points on moduli space where monopoles or dyons become massless are local 
minima of the free energy. There is no evidence for any other local minima. Assuming so, 
this means that at sufficiently low temperatures there are two distinct equilibrium states, 
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related by a spontaneously broken discrete i?-symmetry. In contrast, at sufficiently high 
temperatures there is a unique equilibrium state and the /^-symmetry is unbroken. Hence, 
the theory undergoes a thermal phase transition. The transition temperature must be a pure 
number times the strong scale A. 

Our analysis mirrors that of Ref. [3], but also extends it in several important areas in 
order to fix the misunderstanding of the free energy. The problematic interpretation suggested 
by Figure 1(a) arises from a sign error and a mistreatment of zero- frequency modes. We 
determine the sign of a crucial next-to-leading term in the momentum expansion of the low 
energy effective theory using general arguments based on analyticity constraints satisfied by 
scattering amplitudes in any UV-complete theory [4]. An S'-duality transformation is then 
used to relate results in different regimes of moduli space [5]. When a large hierarchy separates 
the temperature from smaller momentum scales of interest, we construct appropriate three- 
dimensional effective theories which implement the Wilsonian procedure of integrating out 
short distance fluctuations in order to generate effective descriptions for the long distance 
degrees of freedom [6, 7, 8, 9]. 

The remainder of this paper is organized as follows. In Sec. 2 we review relevant facts 
about 577(2) M =2 gauge theory at zero temperature, and discuss the formulation of four- 
dimensional low energy effective theories that will be useful for studying the low temperature 
regime. In Sec. 3, we discuss the high temperature limit and the unique equilibrium state 
that realizes all /^-symmetries. Portions of this analysis involving the construction of the 
appropriate high temperature three-dimensional effective theory are relegated to Appendix A. 
Low temperature thermal effects on moduli space are the subject of Sec. 4. We analyze the 
thermal effective potential for the scalar field that is related to the local coordinate on moduli 
space. The generalization of our analysis to so-called N = 2* theory, obtained by adding 
a single flavor of massive adjoint hypermultiplet to N = 2 super- Yang-Mills, is discussed in 
Sec. 5. Finally, in Sec. 6 we summarize our findings and discuss some open questions. 

2. Review of SU(2) Af=2 gauge theory 

We consider four dimensional N = 2 supersymmetric pure Yang-Mills theory with gauge group 
SU(2). It is renormalizable and asymptotically free. Consequently, the dimensionless running 
coupling transmutates into a renormalization group invariant energy scale A. This theory 
describes the interactions of an N = 2 vector multiplet A. In terms of N = 1 superfields, 
the superfield A consists of a scalar-valued adjoint representation chiral multiplet $ and a 
spinor- valued chiral field strength W a . On-shell, the component fields in $ are a complex 
adjoint scalar (p and an adjoint Weyl fermion ip a . The field strength W a contains an adjoint 
Weyl fermion A Q and a gauge field A^. We shall always work in Euclidean space unless noted 
otherwise, with an action 1 





To obtain a Minkowski space Lagrange density, £(Mmk.)i from the Euclidean version, one performs the 
rotation x° = — ix% and identifies £(Mink.) = — £(Euci.)- 
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Using J\f = 1 superspace notation, the Lagrange density is given by 

-g 2 C= (J d 2 6 \ti{W a W a ) + H.c.^ + 2 j d 2 9d 2 9tr(&e [2V '- ] <S>), (2.2) 

where the field strength W a and vector superfield V are related by 2 

W a = -lD & D«(e- 2V D a e 2V ). (2.3) 

We employ a matrix notation where all fields are Lie algebra-valued (so, for example, W a = 
W®T a with repeated group indices summed over a = 1, . . . , dim G). We take G = SU (N c ) and 
will specialize to N c = 2 momentarily. The fundamental representation Lie algebra generators 
T a are traceless Hermitian N c x N c matrices satisfying [T a ,T b ] = if abc T c , normalized such 
that ti(T a T b ) = S ab /2. The structure constants are real and totally antisymmetric. The 
integrands of the superspace integrals are manifestly invariant under gauge transformations 
of the form e 2V — > e _iA e 2 ^e iA , where A is a fundamental representation chiral superfield. 
Under gauge transformations, the fields W a , 3?, and e' 2,/ ''l < 3? all transform via conjugation by 
the same group element. 

Starting from Eq. (2.2), it is straightforward to show that the Lagrange density in terms 
of on-shell component fields is given by 

g 2 C = 2tr{|F^F^ + iXa^D^X + fya&Drf + (D ^ D ^ 

1 . (2.4) 

- iV2[X, V# f - iV2[X, V# + |[0t, 0] 2 | . 

The covariant derivative = <9 M + i[A^,-], and the field strength F^ u = d^A u — d y A^ + 
z[>1m,j4 v ]. For a Weyl fermion, Aq, = (X a )^ where a, a = 1,2. Our spinor conventions follow 
those of Ref. [10] except that the metric is 5^ and the matrix cr^ is anti-Hermitian. The 
matrices (a^) a a form a basis for 2x2 complex matrices: cr^ = i(~Q _?), and are the 
standard Pauli matrices. The e tensor is used to raise spinor indices to obtain (<7g)^ = 
e /3a e l3a (a^) a a- Numerically, = and = —cr^, although the index structures are 
distinct. Note that [X,tp] = X a ip a - ip a X a and [X,ip] = X^ - tpdX'- 

The Lagrange density (2.4) is invariant under a global SU(2)r x U(1)r i?-symmetry. The 
fermions (^) transform as a doublet under SU(2)r while and <f> transform as singlets. 3 
The U(1)r factor is an ordinary M = 1 ^-symmetry under which <1> has charge 2 and W a has 
charge 1. Quantum mechanically, the U(1)r is anomalous and only a Z4tv c subgroup survives. 4 
For N c = 2, the true global i?-symmetry is thus (SU(2)ji x Zg)/Z2, where the division by Z2 



2 We use the superspace conventions of Ref. [10]. The factor of 2 in front of the Kahler potential ensures 
that the kinetic terms of the two Weyl fermions have the same normalization. 

3 Since A and ip belong to different JV = 1 multiplets, one may check the consistency of Af = 2 supersymmetry 
in Eq. (2.4) by testing the invariance of the Lagrange density under ( ^ ) — > ( jf) • This discrete transformation 
corresponds to the ( _J 0) element of SU(2)r, which generates a Z4 subgroup. 

4 This anomalous i?-symmetry is the reason why no topological charge density (and associated theta angle) 
appears in the Lagrange density. Since $ = <f> + ^/29tp + ■ • ■ , and W a = — i\ a + ■ ■ ■ , it follows that both 
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is a reminder not to double count the (— 1) symmetry (with F fermion number) present in 
both the center of SU(2)r and Z§. 

It will be useful for later purposes to mention another massless representation of J\f = 2 
supersymmetry, the hypermultiplet Ti. In terms of M =1 superfields, % consists of two 
scalar- valued chiral multiplets Q and Q' that transform under conjugate representations of 
the gauge group. On-shell, Q contains a complex scalar q and a Weyl fermion ip q . Similarly, 
Q' contains a complex scalar q' and a Weyl fermion t/y . The scalars ( q \ ) transform as an 
SU(2)r doublet while ip q and ?/y transform as singlets. Both Q and Q' have -R-charge 0. 

Vacua in this theory may be described classically by the requirements that = A = 
ij) = 0, (j) is covariantly constant, and \(j>, <jr\ = 0. If a diagonalizing gauge transformation 
is made to write (j) = aa^/2 for some arbitrary complex number a, then (f> automatically 
commutes with its Hermitian conjugate. The two eigenvalues of 4> are a and —a. Since we 
are free to permute them, <f> = —aa 3 /2 also describes the same vacuum. This permutation 
freedom is part of the residual gauge invariance (specifically, conjugation by ( _? jjj))- 

A translation invariant and gauge invariant "order parameter" parameterizing the space 
of vacua is u = (tr(0 2 )). Since u is a complex number, the classical space of vacua is the 
complex u-plane. In the quantum theory <f> (and its eigenvalues ±a) is a fluctuating field. 
The space of vacua may be changed by quantum effects, but it can never be entirely lifted. 
One reason is that it is impossible to generate an effective superpotential (and therefore no 
squares of F-terms in the scalar potential) invariant under N = 2 supersymmetry without also 
including at least one light hypermultiplet. There are no hypermultiplets at weak coupling. 
For | (a) | 3> A, asymptotic freedom ensures that the theory is weakly- coupled. When quantum 
fluctuations of the eigenvalues of <\> are small compared to their vacuum expectation value, 
the fractional difference ((a — (a)) 2 ) /{a} 2 <C 1 and hence u « (a) 2 /2. 

For any vacuum satisfying \u\ 3> A 2 , a weak-coupling mean field analysis is reliable and 
shows that the Higgs mechanism reduces the gauge group from SU(2) to U(l). There are 
massive W bosons charged under the U(l) photon. The W bosons have masses proportional 
to the expectation value of a, 

M w = V2\(a)\. (2.5) 

J\f = 2 supersymmetry requires the W bosons to belong to Abelian vector multiplets, and 
other components in the multiplet must have the same mass. The dynamics of the resulting 
theory at momenta much less than Myy is both Abelian and N = 2 supersymmetric. The Zs 
i?-symmetry is spontaneously broken since u has i?-charge 4. The unbroken subgroup Z4 acts 
trivially on u, while the coset e 27rJ//8 Z4 acts as u — > —u. 

In Ref. [1], it was shown that when quantum effects are taken into account the space 
of vacua, or moduli space, is precisely the u-plane but with three singular points: one at 

Weyl fermions have 7?-charge 1. If one repackages these Weyl spinors as a single Dirac spinor $ = ( ) , then 
U(1)r acts as a continuous chiral transformation, 'J' — > e'" 75 ^, in a basis where 75 = diag(l, 1,-1,-1). Under 
this field redefinition, the fermion measure in the functional integral acquires a nontrivial Jacobian involving 
the exponential of the topological charge. Therefore, an appropriate choice of w allows one to cancel any 
dependence of the theory on the theta angle. 
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infinity and two at finite values ±uo- One may choose to renormalize the operator tr(c/> 2 ) 
so that uq = A 2 . The existence of a continuous set of vacua implies that (tr (^> 2 (rc))) may 
have arbitrarily long wavelength fluctuations. Such configurations can have arbitrarily small 
spatial gradients with negligible cost in energy, implying that there are massless states in the 
spectrum of the Hamiltonian. These states comprise an M = 2 Abelian vector multiplet. The 
singularity at infinity is a consequence of asymptotic freedom. The singularities at ±uo are 
interpreted as vacua in which extra massless states appear in the spectrum. These massless 
states have spins and 5, and constitute an J\f = 2 Abelian hypermultiplet. Since there are 
no elementary hypermultiplets in the theory, these extra particles are solitonic excitations. 
Massless non- Abelian gluons never appear for any choice of u. That is, there is no vacuum 
corresponding to an infrared fixed point with conformal invariance. Every choice of u (even 
u = 0) corresponds to a theory in which the long distance dynamics is Abelian, possibly with 
extra massless excitations. The theory is always in a Coulomb phase. 

To discuss the particle spectrum at a given u, it is helpful to construct an effective 
theory describing the dynamics in such a vacuum at arbitrarily low momentum. At a generic 
point in moduli space, the massless fields comprise a U(l) N = 2 vector multiplet which is 
simply the neutral component A 3 = ( < 1 )3 , W„) of the gauge triplet. 5 The complex scalar field 
cj) 3 = & 3 \ g= g =0 is identical to the eigenvalue field a when <p is diagonal. Abusing notation 
(but following Ref. [1]), we shall henceforth refer to A 3 as A, its jV=l scalar- valued chiral 
multiplet $ 3 as A, and its field strength W 3 as W a . The low energy effective theory possesses 
N = 2 supersymmetry, and this is made manifest by constructing a Lagrange density directly 
in A^= 2 superspace, 



The prepotential J~(A) is a mass dimension two holomorphic function of A. The function 
JC(A,A) is dimensionless and non-holomorphic in A. e The terms in Eq. (2.6) are organized 
as an expansion in the 'order in derivatives' n, explained in Ref. [5]. The number n is defined 
such that A has n = and the supercovariant derivative has n = 1/2. Ordinary spacetime 
derivatives have n = 1 since they arise from anticommutators of supercovariant derivatives. 
From the structure of a supercovariant derivative, it immediately follows that Grassmann- 
valued superspace coordinates 9f,6 m have n = —1/2. Gauge fields and scalars have n = 
and fermions have n = 1/2. Based on this counting scheme, the chiral superspace integral has 
n = 2 and the full superspace integral has n = 4. Therefore, knowledge of the prepotential 
completely determines terms in the effective Lagrange density with up to two spacetime 
derivatives and at most four fermions. Note that Eq. (2.6) is unchanged by a linear shift of T 

J This is a direct consequence of the Higgs mechanism for large By analytic continuation in the u-planc 
(avoiding possible singularities or cuts), it must also be true even for jitj ~ A 2 where the dynamics is strongly- 
coupled. In fact, if this were not true, then the photon would have to obtain a mass through some type of 
Higgs mechanism. This cannot happen because there are no charged scalars in the 17(1) vector multiplet [11]. 

6 Henceforth, A = J$ . 
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or the addition of a holomorphic function (and its Hermitian conjugate) to fC. When T and 
K. are determined through matching calculations at the momentum scale Mjy, then Eq. (2.6) 
will correctly reproduce gauge invariant correlators of the light fields for distances 3> . 
The effective Lagrange density (2.6) in M = 1 superspace notation is [1, 5] 



where 



and 



r n=2 



Itt 



Im 



cff 



/ 



£2r 2 + CfF 4 + 0(n>6), 



d 2 9 l^" {A) W a W a + / d 2 9d 2 9 T'{A) A 



(2.7) 



(2- 



/~n=4 



'cff' 



(2.9) 



d 2 9 d 2 9 [lC AA {A, A) (D a D a A)(D & D & A) + 2(D & D a A){D a D & A) 
+ 4(D a W a )(D 6l W^) 

- A{D^W^){D [a W p) ) - 2D a D a (W W ) 

- A(D {6i W^)(D^W^) - 2D & D & (WpW^) 

- 2IC AAA (A, A)W a W a D?DpA - 2IC AAA (A, A)WaW & D^A 

+ KaaAa(A A) [-8(W a D a A)(WaD«A) + AW a W a W & W & ] } ■ 

In expression (2.9), subscripts on the non-holomorphic function K, indicate derivatives with 
respect to the indicated arguments. Expression (2.9) is unique up to terms proportional to 
D a W a — DaW a which vanish because W a satisfies the Bianchi identity. 

A gauge invariant description of moduli space is the u-plane, a one-complex dimensional 
manifold. One may promote the coordinate u to a field u(x), valued in the complex numbers. 
The dynamics of arbitrarily long wavelength fluctuations of u(x) around a constant value 
is described by a sigma model action of the form S e g = J d 4 x [y(u, u) d v ud v u + •••]. One 
can regard the coefficient of the two derivative term as a metric on moduli space. The line 
element is written as ds 2 = j(u, u) du du. 

The metric on moduli space is easily determined for \u\ S> A 2 . In this regime the effective 
theory is formulated in terms of the scalar field a(x), and its interactions are weakly coupled 
due to asymptotic freedom. The translationally invariant vacuum expectation value (a) is 
mapped to u by the approximate formula u ~ (a) 2 /2. Therefore, the asymptotic region of 
moduli space may be parametrized by the local coordinate a. 7 The induced metric in field 
configuration space is obtained from the Kahler potential K(A,A) = ^Im(J-'(A)A). The 
full superspace integral of the Kahler potential yields 



C 



cff 



-f(a,a)d fl ad ll a + 



with the Kahler metric 



7(0, a) = d a d a K\ g= 



(2.10) 



(2.11) 



7 Here a is understood to mean (a). Future usage should be clear from context. 
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However, the effective theory is more than just a sigma model; it is also an Abelian gauge 
theory. Define a holomorphic gauge coupling function t(A) = J 7 " (A). The half superspace 
integral of ^t(A) W a W a yields 

C eS = \g^(a, a)F^F^ u + -U, 6 cS (a, a)F^F^ u + ■■■ , (2.12) 

where the inverse effective gauge coupling and the effective theta angle are given by 

g~$(a,a) = ilmr(a), O cS (a,a) = 2nRe r(a). (2.13) 

Note that M = 2 supersymmetry requires that the inverse gauge coupling and the Kahler 
metric 7 coincide. 

The prepotential fixes all the effective couplings in C^ 2 . It is straightforward to express 
expression (2.8) in terms of off-shell component fields, 

Cf 2 = 9~£(a, «){|^a| 2 + + (fyo&Drf + f Aa^A + H.c) - |F| 2 - ±D 2 

+ [r(a, a) (^'i 7 * + jA 2 F - -^XipD - ^Xa^^F^) + H.c] (2.14) 
- [i?(o, a)^AV + H.c] } + 35^ ^ C ff(a, a) F^ U F^ U . 

In expression (2.14), each fermion bilinear is shorthand for a spinor contraction (e.g., ■0 2 = 
^"Va), the Abelian field strength F^ = d^A y - dyA^ and its Hodge dual W v = \^ vafi F afi 
(with e 0123 = 1), and 

Dp = r(o, Ot) djjfl j (2.15a) 

T(a,a) = 7(0, a) -1 9 a 7(a,a) , (2.15b) 
JR(a, a) = <9 a r(a, a) + T(a, a) 2 

= 7 (a,a)- 1 a 2 7 (a,a). (2.15c) 

The on-shell form of O^ 2 may be obtained by solving the equations of motion for the auxiliary 
fields, but we will not need that result. 8 

The prepotential may be determined as follows. The one-loop beta function for the 
running gauge coupling of SU(2) M = 2 gauge theory is 



dg 2 1 4 



l + 0(e- 87r2 / 9 ' 2 ) . (2.16) 



Higher order loop corrections vanish due to supersymmetry, but we have included the form 
of nonperturbative one-instanton corrections. Integrating Eq. (2.16) yields 

g- 2 ^) = _L hV/A 2 ) + const. + 0(A 4 //) , (2.17) 

4"7T Z 



8 The D-term equation is D — — ^T(a, a)i\ip + H.c. and the F-term equation is F — §r(a, a)ip 2 + 
|r(a, a)*A 2 . Substituting these into expression (2.14) produces four-fermion operators only. We do not need 
the explicit result since we will carry out perturbative calculations using Feynman rules for off-shell fields. 
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where A is the conventional definition of the strong scale. This may be matched to the effective 
gauge coupling of the low energy effective theory at the W mass scale. That is, at fj, = My/ 
one has g 2 s (a,a) = g 2 (Mw)- For asymptotically large \u\, the mass formula Eq. (2.5) implies 



The line element on moduli space is ds 2 = jrlm r{a)dada, where the metric is given 
explicitly by ^lmr(a) ~ |^[ln(|a| 2 /A 2 ) + 3]. The metric is single-valued and positive for 
\a\ S> A. It diverges as \a\/A — > oo which means that the effective gauge coupling becomes 
arbitrarily small. This is just a restatement of asymptotic freedom. For smaller values of \a\ 
there is a difficulty: the metric is negative. In fact, r(a) is holomorphic so Im r(a) must be 
harmonic, and since it is not constant, it must be unbounded from below. The metric fails 
to be positive-definite, or equivalently, g c f[ fails to be real. One is forced to concede that a is 
valid as a local coordinate only asymptotically far out on moduli space. A key observation 
of Ref. [1] is that the form of the metric, as well as its positivity, can be maintained if an 
additional coordinate on moduli space is introduced: ao = dJ-(a) / 'da. Then the line element 
on moduli space may be expressed as ds 2 = ^(dap da — dadar,) which exhibits a symmetry 
under ( a g ) ->• ( -l D ) . This allows one to also use an as a local coordinate on moduli space, 
with a different harmonic function serving as the metric. The region of the u-plane in which 
ap is a good coordinate will be discussed shortly. 

The metric on moduli space is preserved by real linear fractional transformations of the 
coordinates ( a <?). Define a holomorphic vector field a(u) = ( a ^^) over the u-plane. The 
induced metric is given by ds 2 = —^e mn ^^^^dudu, and is preserved by monodromies 
M G SX(2, Z) that act on the vector field as a — > Ma. 9 

The presence of monodromy and M = 2 supersymmetry naturally imply a type of electric- 
magnetic duality. This fact is uncovered by asking the following question: if monodromies 
rotate a into ao, but a belongs to a vector multiplet, then what effect does an SL(2,Z) 
transformation have on A^l In particular, consider the monodromy 



which fully rotates a into — a/> In Ref. [1], it was demonstrated that S acts on the gauge 
fields as a Fourier transformation in field configuration space from to a dual gauge field 

9 The invariance of the mass formula for BPS-saturated states under the action of monodromies implies that 
M must be integer-valued, and that no constant vector can be added to the linear transformation a — >• Ma. 
The actual monodromy group turns out to be T(2) C SL(2,Z) which consists of matrices congruent to the 
identity matrix (modulo 2, taken element-wise) [1]. 




(2.18) 



This leading log is reproduced by a prepotential 




(2.19) 




(2.20) 
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Ad^l- The form of the action remains the same in the new variables except that the effective 
gauge coupling inverts (as a consequence of integrating a Gaussian functional). In summary, 
SL(2, Z) acts linearly on the J\f = 1 chiral multiplets A and Ad = dJ-(A)/dA, and by electric- 
magnetic duality on the M = 1 chiral field strengths W a and Wn a . Just as -C^f 2 contains 
two derivative terms for the effective theory asymptotically far out on moduli space, one can 
define an analogous Lagrange density for the S'-dual effective theory, 



rn=2 _ l T 
L -D,eS — ~ 4^ 1111 



J d 2 6\F' D {A D )W£W Da + j d 2 9d 2 6r D (A D )A D 



(2.21) 



This description is useful precisely where the original gauge coupling blows up. Thus, the 
S-dual effective theory is valid near the massless monopole point in moduli space. The dual 
prepotential Td may be related to the original prepotential via a Legendre transform [5]. In 
practice, one may obtain .Fd(Ad) using a method similar to the one used to obtain J-{A). 

Four derivative terms in the S'-dual effective theory are given by the Lagrange density 
^zTcff- The explicit form for C D = ^ S is identical in structure to expression (2.9) except that dual 
chiral superfields Ad appear in place of elementary ones, and the dynamics is determined by 
a non-holomorphic function tCp. A priori, there is no relation between K, and ICo- However, 
in Ref. [5] it was proved directly in M = 2 superspace that the non-holomorphic function tC 
is a modular function with respect to SL(2, Z). In particular, under an S transformation 
(which amounts to a Fourier transformation), 

K D (-A D ,-A D ) = 1C(A,A). (2.22) 

This relation will be important later. 

The Af = 2 theory has a BPS bound, M > y/2\Z\, where Z is the central charge in the 
extended supersymmetry algebra. The lightest states saturate this bound. A state with 
electric charge n e (defined by its coupling to the photon A^) and magnetic charge fi m has 
mass 

M = V2\Z\, (2.23) 

where Z(u) = a{u)n e + au(it)n m . Since Z determines particle masses, it is renormalization 
group invariant. 

The singularities in the finite ti-plane arise from massive J\f = 2 Abelian hypermultiplets 
that become exactly massless at u = ±uo [1]. At uq, a magnetic monopole with charges 
(n m ,n e ) = (1,0) becomes massless. At — uo, a dyon with charges (n m ,n e ) = (1, —1) becomes 
massless. 10 Let us study the monopole. According to Eq. (2.23), the monopole has a mass 
proportional to the coordinate od, 

M m = V2\a D \. (2.24) 

10 In Ref. [1] these charge assignments for the extra massless multiplets were shown to pass several consistency 
checks. For instance, there is a monodromy around each singularity in the u-plane and the set of monodromies 
should furnish a representation of the fundamental group of the u-plane (with singularities deleted) in SL(2, Z). 
If the monodromies around «o and —Mo are calculated assuming that the hypermultiplets becoming massless 
are a (1,0) monopole and (1,-1) dyon, respectively, then the monodromies indeed generate a subgroup of 
51/(2, Z). As another check, the triangle inequality and conservation of energy ensure that when the ratio an I a, 
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This mass vanishes at uq only if cld{uq) = 0. In the vicinity of the point uq, the low energy 
effective theory must describe a U(l) M = 2 gauge theory coupled to a massive hypermultiplet. 
This theory is essentially an J\f = 2 version of QED with the light monopoles playing the 
role of "electrons." It is infrared free and the renormalization of the dual effective gauge 
coupling is due primarily to one-loop photon self-energy diagrams where the fields of the light 
hypermultiplet run around the loop. Consequently, the renormalization group implies that 
each decade of momentum, from a fixed UV cutoff down to the monopole mass, contributes 
the same amount to the inverse coupling. Since the monopole mass vanishes as u — > uq, it 
follows that the dual effective gauge coupling gp vanishes as u — > uq. Hence, is a good 
coordinate on moduli space in a neighborhood of the point uq. 

To determine the dual prepotential, consider the beta function of N = 2 QED with a 
single hypermultiplet, 

Integrating Eq. (2.25) yields g] D 2 (M m ) = ^ ln(A 2 /Af^) + const. This result holds for u close 
to uq. The inverse gauge coupling is related to the imaginary part of a holomorphic function 
t d {A d ) = T'^Ad). Plugging in Eq. (2.24) yields 

ImT D (a D )K-±ln(^pj. (2.26) 
This leading log is reproduced by a dual prepotential 

fa 2 
,2 ,„/ a D 



F D {a D )^-—a' D ]n[-^\. (2.27) 

An explicit expression for the four derivative terms in C^ 4 , including a formula for the 
non-holomorphic function JC, is discussed in Sec. 4. 

Because there are i?-symmetry transformations which act on the entire vacuum manifold 
as u — > —u, the behavior of the free energy density near the massless monopole and dyon 
points is identical. Consequently, it is unnecessary to write down a low energy effective theory 
incorporating light dyons. Our low temperature analysis relies only on the weakly-coupled 
effective descriptions near u = oo and u = uq. 



3. High temperature behavior 

SU(2) M = 2 gauge theory at asymptotically high temperatures, 

T>c/T>c/ 2 r>A, (3.1) 

is not real, the monopole and dyon cannot decay because n m and n e are relatively prime. Asymptotically 
far out on moduli space, there exist field configurations with magnetic charge — these are the semiclassical 
monopoles. Moving in from infinity along the real u-axis toward ±uo, one never crosses a curve on which an /a 
becomes real. This implies that whatever stable BPS-saturated states exist at infinity must also appear in the 
strongly-coupled region. 
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is weakly coupled, g -C 1, on length scales small compared to 1/(<? 2 T). [Here and henceforth, 
g = g{T) stands for the running gauge coupling evaluated at the scale T.\ Given the hierarchy 
of scales (3.1), one may compute the dependence of the free energy density F/V on the 
(translationally invariant) thermal expectation value of the complex scalar field (f>, using 
effective field theory techniques and perturbation theory. In other words, one may compute 
the thermal effective potential for <p. Minimizing F/V with respect to (<j)) determines the 
number and location of equilibrium states, which in turn determine the realization of the 
discrete i?-symmetry. 

High temperature perturbation theory, at ((f)) = 0, shows that the non-Abelian N = 2 
plasma has a positive Debye mass, 



Since the curvature of the thermal effective potential for (f>, at ((f)) = 0, equals the thermal 
mass m?, the positive value (3.3) indicates that (p = is a local minimum of the free energy. 
To demonstrate that this is, in fact, the global minimum, one must evaluate the effective 
potential for (p arbitrarily far away from <p = 0. This we do in Appendix A. The result 
is unsurprising: for asymptotically high temperatures there is a unique equilibrium state at 
((p) = 0. The free energy density in this equilibrium state has the asymptotic expansion 



The leading term is the ideal gas blackbody contribution. The 0(g 2 ) term comes from 
two-loop contributions at the momentum scale of T, while the 0(g 3 ) term arises from zero 
frequency contributions on the scale of gT. 

The positive thermal mass (squared) (3.3) and the unique minimum of the thermal ef- 
fective potential imply that the discrete i?-symmetry, which is spontaneously broken at zero 
temperature, is restored for sufficiently high temperature, T ^> A. 11 To argue this formally, 
recall that the gauge invariant order parameter involves the operator tr ((f) 2 )- To probe sponta- 
neous symmetry breaking one may consider the correlator D(x—y) = (tr ((p(x)^) 2 tr ((p(y) 2 )), 
where (...) represents an expectation in a .R-symmetry invariant thermal equilibrium state 
(which might be a statistical mixture of two noninvariant pure states). A non- vanishing large 
distance limit, D(x—y) -/->0 as \x— y\ — > oo, would indicate breakdown of cluster decompo- 
sition and consequent spontaneous symmetry breaking. However, the positive thermal mass 
(3.3) implies that scalar correlators, evaluated at the ((f)) = global minimum of the effective 
potential, fall exponentially fast at distances large compared to rri^ . Consequently, D(x — y) 
approaches zero at large distances and does not contain a disconnected part, signaling unbro- 
ken R symmetry. 

11 A similar conclusion was reached in Ref. [14]. However, this analysis did not include the effects of inter- 
actions. 





and that the field cp also develops an effective thermal mass 




(3.3) 




(3.4) 
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4. Low temperature effective theory 

The spectrum and low energy dynamics of SU{2) N = 2 gauge theory varies drastically over 
moduli space. Far out on moduli space, near the asymptotically free vacuum, charged particle 
masses become arbitrarily large and the low energy dynamics contains only an Abelian mass- 
less vector multiplet. Near the strongly-coupled vacua there are, in addition to the massless 
vector multiplet, charged hypermultiplets consisting of light magnetic monopoles or dyons. 
Whether or not this extra matter affects the low energy dynamics is determined by the ratio 
of the monopole (or dyon) mass to the temperature. We will focus attention on three distinct 
regions of the free energy surface where weakly-coupled effective theories may be constructed. 
Using these effective theories, our goal is to compute the free energy density as a function of 
the moduli space coordinate u. 

4.1 Near the singularity at infinity 

Consider the asymptotically free region of moduli space. Near u = oo, the spectrum of the 
theory includes very heavy BPS states satisfying M m S> Mw 3> A. Suppose the temperature 
T <C M\y. In this regime, it is valid to use an effective description of the low energy physics 
in terms of the scalar field a(x) discussed in Sec. 2. The appropriate effective theory contains 
only the U(l) vector multiplet A and is described by the Lagrangian C e g given in Eq. (2.7). 

Although it is possible to construct an effective theory incorporating the additional 
massive charged vector multiplets (see, for instance, Ref. [1]), the W bosons have mass 
M\y = y/2\a\ which, by assumption, is large compared to T. The Boltzmann weight of 
these W bosons exponentially suppress their contribution to the free energy density relative 
to the contributions of the massless particles. Since the W bosons and their superpartners 
form a dilute gas, it is straightforward to obtain 12 

(F/y) chargcd «-16T 4 (^) 3/ V^/^. (4.1) 

Since the pressure is just minus the free energy density (when all chemical potentials vanish), 
the negative sign in this result shows that the dilute gas of heavy particles exerts a positive 
pressure, as it must. Eq. (4.1) also indicates that the pressure decreases as Mw/T grows large. 
This may be achieved, at fixed temperature, by moving toward large \u\ in the free energy 

12 Each charged degree of freedom in the plasma generates a contribution to the free energy density equal 
to /± = ±|Tj] ngz / jl^V ln[(w ? t) 2 +p 2 + M w ], where u+ = 2nnT for bosons, and ui~ = (2n+l)nT for 
fermions. The mass Mw is the same for all particles belonging to a given multiplet. Note that Mw cuts 
off the infrared divergence in the spatial momentum integral of the bosonic zero-frequency contribution. The 
temperature-dependent part of the free energy contribution f± is finite, and is a standard integral in statistical 
mechanics, f± — (T-indep.) ± J Q °° dx x 2 ln(l =p e ~V x2+M w/ T2s j . In a charged vector multiplet, the on-shell 
bosonic fields (gauge field and complex scalar) account for four real degrees of freedom. By supersymmetry 
the on-shell fermionic fields (two Weyl fermions) must also account for four real degrees of freedom. There 
are two such charged multiplets for SU(2) broken to U(l). Thus, (-F/V) c harged = 8 (/+ + /-)• Evaluating the 
leading contribution from the saddle-point at x = immediately gives the result (4.1). 
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surface. Also, note that the contributions of monopoles and dyons to the free energy will have 
even greater exponential suppression since these excitations, for large \u\, are much heavier 
than W bosons. This follows from the mass ratio of electrically and magnetically charged 
BPS-saturated states near infinity on moduli space, M m /M\y = \ajj\/\a\ ~ | ln(u/A 2 )|. Thus, 
as noted in Ref. [3], one may ignore massive charged fields altogether and focus on the 
interactions of the massless neutral fields. 

To compute the effective potential V e s for the scalar field a, one may first integrate 
out thermal fluctuations of the vector multiplet fields in perturbation theory. As usual, one 
considers the functional integral representation for the partition function, with a constant 
source coupled to the scalar field, and performs a saddle point expansion. Let ao represent the 
translationally invariant expectation value for the scalar field in the presence of the source. 
By construction, ao is a solution to the Euler-Lagrange equations for a fixed value of the 
source. Let 

a(x) = a + a(x). (4.2) 

Substituting this decomposition into the Lagrange density £™f? 2 [gi yen by expression (2.14)] 
and using the prepotential (2.19) leads to 

r n=2 _ AO) , Al) , A 2 ) j (AO) 

The coupling constant g 2 . = 4-7r 2 /(ln |ao/A| 2 + 3) is small in the asymptotic regime under 
consideration. The perturbative expansion for V e R will be a series controlled by g$. The free 
Lagrange density is 

Cf cc = \d»~a\ 2 + \Fl v + {l^d^ + §A^A + H.c.) - |F| 2 - \D\ (4.4) 
The interactions cubic and quartic in fluctuations are given by 

^^ = (A + ?!) £ w + i(i-f)a.+(f<t-+^). 

g \ ao ao 7 4 \ ao ao / v ao / 

4^ 2 (2 )_ I fa 2 5* 2 \ (o) if a 2 a* 2 \ ~ / a 1 
~gj C ^ " "2 lag + 4) Ci ™ ~ 8 la| " ^) F ^ + \~4 ° fcrmi + 4a| A ^ + Hx ' 

(4.5b) 

where 

Ofcrmi = -fOKV 5 + + \^ 2 F* + \\ 2 F - -^\i[)D - ^AOi^ (4.6) 

is a dimension five operator composed of fermion bilinears. To obtain Eqs. (4.4) and (4.5), all 



fluctuating fields have been rescaled by a common factor of g$. In general, the contribution 

(p) 

£- m { to the interaction Lagrange density involves 2 + p factors of fluctuating fields. Each such 
term has an overall factor of (?Q +2 /|ao| p multiplying operators of dimension 4 + p. 13 



13 The reason for p > 0, contains an overall factor of g^ +2 (instead of just g^) is due to the fact that in 
the original expression (2.14) for £™Jf 2 , the Kahler connection T(a, a) and R(a, a) contain an inverse power of 
the Kahler metric j(a, a) which (because 7 coincides with g~g) cancels the overall factor of g~g. 
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A schematic expression for the effective potential is 

r 2 . / ^\1PI / / o\ . /,v„\lPI / /on m ™ . ,,v,a1PI 



The first term in Eq. (4.7) is the blackbody contribution from a massless M = 2 vector mul- 
tiplet. 14 In the remaining terms, we have omitted spacetime integrals and combinatorial 
coefficients. A generic term of the form (n^=i ^-int^)o PI represents the expectation value of 
m spacetime integrals of the interaction Lagrange densities in the unit-normalized Gaussian 
measure. If this term is expressed diagrammatically, then each diagram will have m vertices 
(representing insertions of the C^l ), joined by propagators arising from C^ e . Only the one- 
particle irreducible (IPI) portion of these correlators contributes to the effective potential. 

The effective potential admits a double series expansion in the dimensionless coupling go 
and the ratio of scales T/\ao\. To this end, the expansion in Eq. (4.7) has been organized by 
operator dimension, with the dimension 5, 6, 7, and 8 terms shown explicitly. The Gaussian 
measure is invariant under independent U(l) phase rotations for a, ip, A, and F, and Z2 parity 
transformations for and D. These symmetries immediately imply that the dimension 5 
and 7 terms in Eq. (4.7) (which come with odd powers of go), and the term (jC^|)J vanish 
identically. 

From ((^int) 2 )^ 1 there are eight basic diagrams that must be considered. These "bas- 
ketball diagrams" are shown in Figure 2. 15 All of these diagrams vanish because they reduce 
to the sum-integrals 16 



tp,± Vq,± P Q 

or 



1 = 0. (4.8b) 



± 



Expression (4.8a) vanishes by Euclidean time reflection and spatial parity invariance. To 
justify this one may choose to regulate the theory by dimensional continuation, which pre- 
serves spacetime symmetries. [Or one may ignore the issue of regulation, since that is part 



14 Note that there are no background-dependent terms in C^ ee . Consequently, the one- loop contribution to 
V e fi (involving the logarithm of a functional determinant) simply gives the blackbody result. 

15 Although auxiliary field propagators are momentum-independent, the diagrams containing F and D prop- 
agators have two independent loop momenta. This is clear if the diagrams are constructed using position-space 
Feynman rules. Diagrams 2c and 2d may be thought of as originating from four-fermion interactions in the 
on-shell formalism. The auxiliary field is a constraint that causes these graphs to "pinch," yielding a graph 
with a figure-eight topology. Such an on-shell diagram has two fermion loops which contribute two overall 
factors of —1; this is matched in the off-shell diagram by one fermion loop and one central auxiliary field 
propagator that each contribute a factor of —1. 

16 Sum-integrals are defined as |^ ± = T ^2 to j ' (^r) 3 ' wnere the sum is over even (+) or odd (— ) integer 
multiples of nT. 
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Figure 2: Two-loop 1PI diagrams that could contribute to the effective potential at 0(g®T Q /\a$\ 2 ). 
Solid lines represent complex scalar fields, dashed lines represent either type of Weyl fermion, wavy 
lines represent Abelian gauge fields, and dotted lines represent auxiliary fields. 



of denning the theory at zero temperature, and focus only on the temperature-dependent 
part of the sum-integral. Each discrete frequency sum may be recast as a pair of contour 
integrals just above the real axis: one temperature-independent and the other temperature- 
dependent. The latter contains the appropriate statistical distribution function which dies 
off exponentially fast in the upper half complex plane. The temperature-dependent piece is 
finite and vanishes by the symmetry arguments.] Expression (4.8b) has a scale-free spatial 
momentum integral that vanishes in dimensional continuation. [Alternatively, in the contour 
integral method, the temperature-dependent integrand is analytic in the upper half plane, 
so closing the contour there produces zero.] Thus, there is no 0((7q T 2 /|ao| 2 ) contribution to 
Vcs/T 4 '. This is consistent with the analysis of Ref. [3]. 

Continuing on to the subsequent terms in the expansion (4.7) for the effective potential, 
one finds that contributions to V^j/T 4 from dimension 8 operators debut at 0(^0 ^ 4 /l a o| 4 )- 
The local dimension 8 piece (£^)o displayed in the Eq. (4.7) has an overall factor of but 

(2) 

(just as for (£- nt )o this piece vanishes due to phase rotation symmetries. 

It was pointed out in Ref. [3] that the first non-vanishing correction to the effective 
potential does not come from the n = 2 terms in the low energy effective action, but rather from 
the higher derivative n = 4 terms. In other words, £"ff 4 and the non-holomorphic function 
K, is more important than C™g 2 and the holomorphic prepotential J-, when it comes to 
understanding the leading dependence of the effective potential on the expectation value of 
the scalar field. To see this one must use the known form of JC [15, 16], 

with c a constant. This is the leading approximation for the non-holomorphic function /C 
when a = A\ e= Q =0 S> A. 17 A couple points regarding the formula (4.9) for K are worth 

17 This term is responsible for producing four derivative terms in the effective action, such as F 4 . (It is 
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Figure 3: Vertices representing leading four-derivative scalar self-interactions. 

mentioning. First, the concise form (4.9) is due to the fact that the low energy gauge group 
is Abelian. 18 Second, both U(1)r and scale transformations (of the form A — > bA) change 
Eq. (4.9) additively by terms that are (anti)holomorphic in the chiral superfield. Such terms 
vanish under the full superspace integration, and therefore have no effect on the dynamics. 

For simplicity, we compute the contribution to the free energy density from the purely 
scalar terms in £ef? 4 - Such terms originate from the piece 

£Sf 4 D-J d 2 6d 2 9IC AA (A,A) [(D a D a A)(D & D A A) + 2(D dt D a A)(D a D«A)} . (4.10) 

We write expression (4.10) in components, perform the expansion (4.2), and rescale fluctuating 
fields by go- 19 The resulting purely scalar terms that are both invariant under phase rotations 
of the fluctuation a, and suppressed by no more than |ao| -4 , are given by 20 



r n=A _ 



^ (|<W 2 + \d 2 a\ 2 ) + 7^4 { (I VI 2 ) 2 + \a\ 2 \d,d u ~a\ 2 + l«| 2 |^ 2 a| 2 



+ [5 (fy&aXfya'X&a*) + H.c] } 



(4.11) 



up to total spacetime derivatives. It will turn out that only the first two terms inside the 
curly braces generate a non-zero contribution to V c g/T 4 of 0(g 4 T 4 /|ao| 4 ). 

The interaction terms (4.11) generate quadratic and quartic interaction vertices, illus- 
trated in Fig. 3, with momentum-dependent vertex factors. Explicitly, the resulting (Eu- 



simpler to discuss perturbative corrections to J d 4 6d 4 6 fC, rather than K, itself, since the former is a Lagrange 
density and thus a direct output of renormalizing the short distance theory.) It has been shown, at least in the 
case of AT = 2 QED, that there is a nonvanishing quantum correction at two-loop order to the supersymmetric 
completion of the F 4 term [17]. In our case, the F 4 term is given by —AC J d 2 9 d 2 9 K-aaaa W 2 W 2 , and the 
corresponding statement about quantum corrections to this operator may be expressed as C = 1 + O(go). The 
loop expansion parameter in the short distance theory is the running coupling g 2 , evaluated at the scale a, 
or in other words, the inverse logarithm (In |a/A| 2 ) -1 . Far out in moduli space, this coupling is small. Such 
two-loop corrections lead to 0((?oT ,8 /|ao| 4 ) contributions to the effective potential which will turn out to be 
sub leading. 

18 The leading form of K, for non- Abelian gauge group SU(2) was originally determined in Ref. [15]. On the 
Coulomb branch of 5(7(2) M = 2 theory, the leading form of K, simplifies considerably, as noted in Ref. [16]. 

19 When chiral superfields occur in denominators, we factor out the lowest (scalar) component and expand 
the remaining function as a power series in 6 and 9. 

20 Contributions from multi-point correlators of terms which are not individually invariant under phase 
rotations of o first contribute to the effective potential at 0(g^T s ' /\ao\ 4 ). This is smaller by two powers of go 
than the contributions which will result from the unperturbed expectation of the [/(l)-invariant terms (4.11). 
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Figure 4: Leading scalar contributions to free energy density up to 0(3qT 8 /|oo| 



clidean space) vertex factors are 

-2\2 



V 2 



^1 ( P 2 r, 



\ a o\ 

m r 

\ a o[ 



(4.12a) 



2(pi ■P2){V3-Pa) + 2(pi -j?4)(P2 -P3) + P1P2 + P2P3 + P3P4 + P4P1 
+ (Pi • P2) 2 + (P2 • Pa) 2 + (pa • P4) 2 + (P4 • Pi) 2 (4.12b) 
+ 2(pi • P2)(pi • P4) + 2(p 2 • P3)(P3 • Pa) + 2(pi • P2XP2 • P3) + 2(pi • £>4)(P3 • P4) ■ 

The momenta are taken along the arrows (which distinguish a* from a). The overall signs 
of both vertex factors are positive because there is a minus sign from the definition of the 
Euclidean path integral weight, a minus sign from the Lagrange density itself, and two factors 
each of i and —i from Fourier transforms of derivatives. 

These quadratic and quartic vertices generate the bubble diagrams shown in Figure 4. 
Diagrams I\ and I2 vanish, since their spatial momentum integrals are scale- free. The last 
diagram, I3, reduces to the square of a nontrivial one-loop sum-integral, 21 

21 o _ _4 / i- z+1 \ 2 



2 typ,+ / Tq,^ 



4 

l«o| 4 



4(p • q) 
p2g2 



15V 



3|a | 4 



(4.13) 



where we have used rotational symmetry to replace fo^l by (1 + gjr[) ^2^2 inside the integrals 
and d is the spacetime dimension (continued infinitesimally away from 4). In Eq. (4.13) the 
sum-integral in parentheses J+ = S ^ evaluates to 7r 2 T 4 /30. The scalar self-interactions 



contribute — (/3V) 1 /3 to the effective potential, so 



27r 4 <7 4 r 4 



V e s(ao) 



scalar-scalar 



rp4 



(4.14) 



675 I a 1 4 

Obviously, the sign of the coefficient c is of the utmost importance for interpreting the effective 
potential — the sign determines whether local equilibrium is directed toward smaller or larger 
values of I czo I - 



21 The last eight terms in the Feynman rule for V% in Eq. (4.12b) contribute terms to the expression for 73 that 
integrate to zero. These eight terms originate from the operators |a| 2 |9 2 a| 2 and a(dfj,d v a)(d^a*)(d v a*) + H.c. 
in expression (4.11). Since these operators have no bearing on the calculation, they are omitted in Eq. (3.6) 
of Ref. [3]. 



- 19 - 



An attempt has been made to derive the non-holomorphic function KL for SU(2) J\f = 2 
gauge theory and fix the value of the overall constant [18]. The method involves technical 
superfield calculations and is difficult to check for errors. 22 Our goal, in the next couple 
of pages, is to find an independent determination of the sign of c, since that is the crucial 
information needed to understand low temperature thermodynamics in this theory. 

A simple and physical method that fixes the sign of the coefficient c is provided by 
studying the forward amplitude for scalar scattering [4]. There are spinless, one-particle 
states in the spectrum of the theory that arise from quantum fluctuations of the field o. Let us 
denote the particle excitation by <p and its antiparticle by (p. Consider the scattering process 
pip —7- ip<p. At center-of- momentum energies far below Myy, one may use the low energy 
effective action to reliably compute the scattering amplitude in a momentum expansion. The 
tree-level diagram for this process is derived entirely from the vertex V2 shown in Figure 3, 
remembering that the Minkowski space vertex gets an extra factor of i relative to its Euclidean 
counterpart. 23 It is worth noting that, from the point of view of the expanded low energy 
effective Lagrange density, the operators appearing in expression (4.11) comprise only a subset 
of the irrelevant interactions. Moreover, they are not even the ones with lowest dimension. 
Nevertheless, the quartic operators in expression (4.11) generate the leading contribution to 
the (pip — > pip scattering process. Contributions from other terms in the effective theory are 
suppressed by additional factors of the dimensionless coupling go- The resulting Lorentz- 
covariant scattering amplitude is 24 



22 In the case of SU (2) N — 4 gauge theory, three independent superfield calculations of the non-holomorphic 
function K, are known, and they completely agree [19, 20, 21]. The calculations for the superconformal theories 
discussed in these works may be readily extended to 5*17(2) M — 2 gauge theory. For example, Eq. (4.9) in 
Ref. [21] represents the one-loop effective action for M — 2 gauge theory with four fundamental hypermultiplets 
and gauge group SU(2) Higgsed to (7(1). Of direct relevance to this work is the sum of the second and third 
terms on the right hand side of Eq. (4.9) — this coincides with the non-holomorphic contribution to the effective 
action for AT = 2 theory without matter. The F 4 term may be extracted by Taylor expanding the functions 
Citfy, tty) and uj{fil , t^) to zeroth order around ^ = = 0. This amounts to focusing on just the leading term 
in a derivative expansion. Since C(0>0) = 1/12 and cj(0, 0) = 0, it follows that c is positive. A positive value 
for c agrees with the conclusion in this work. We thank the JHEP referee for explaining this to us. 

23 For scattering we take all momenta as incoming (rather than following the arrows) and label them as pt, 
i = 1, . . . , 4 starting at the upper left corner and continuing counterclockwise. This amounts to flipping the 
signs of P2 and P4 in Eq. (4.12b), which clearly does nothing to the whole expression. We use a Minkowski 
space with — h++ signature, and define the usual Mandelstam variables s — —{pi +P2) 2 and t = — (p2 +P3) 2 ■ 
The mass-shell condition is p i = 0. 

24 Recall that the LSZ reduction formula relates the Fourier transform of i times time-ordered correlation 
functions to the scattering amplitude M. Thus, —iM is the object to which diagrammatic rules apply. 




(4.15) 



In the forward scattering limit 




(4.16) 
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Figure 5: Analytic structure of the forward amplitude A(s) as a function of (complexified) center- 
of-momentum energy squared. J\f — 2 gauge theory is recovered by sending m gap — > 0. 

Now consider the contour integral 



The analytic structure of the exact forward scattering amplitude in the complex s-plane is 
shown in Figure 5. A(s) must have a branch cut along the positive, real s-axis with a branch 
point that corresponds to the threshold for pair production. Since N = 2 gauge theory has 
excitations with arbitrarily low momentum, there is no mass gap and the branch point sits 
at the origin. Following Ref. [4] , one may modify the theory in the deep IR by giving a small 
regulator mass W/gap to the (x fields. The cut then extends only down to (277ig a p) 2 . There is no 
cut along the negative real axis, since the amplitude M is only symmetric under interchange 
of s and t, not under interchange of s and u. By construction, the integrand „4(s)/s 3 also has 
a pole at the origin. 

The integral / may be evaluated by deforming the contour 7 usefully in one of two ways: 
in a tight circle around the origin yielding the residue at the pole at the origin, or around 
infinity. In the latter scenario the integral along the large circular portion of the contour 
vanishes since SU(2) M = 2 gauge theory is UV-complete and thus the forward amplitude 
(when the theory has a mass gap) grows no faster than s 2 at high energies. This leaves only 
the contour wrapping the cut which measures the integrated discontinuity of A(s) across the 
cut. Since A(s) is real along part of the real s-axis, the Schwarz reflection principle relates 
the discontinuity to the imaginary part of A(s) just above the cut. Consequently, 



On the left-hand side of Eq. (4.18), one may approximate the forward amplitude at weak 
coupling with the tree-level formula. Since the introduction of a gap modifies the mass-shell 




(4.17) 




(4.18) 
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ss gg sg y/yj ' )J. ' y//. sy/ ' gy/ ' gk 

Figure 6: Two-loop diagrams that contribute to the effective potential at O^qT 8 /|a 1 4 ) - Solid lines 
represent complex scalar fields, dashed lines represent either type of Weyl fermion, and wavy lines 
represent Abelian gauge fields. 



condition, the tree- level amplitude is given by Eq. (4.16) with additional terms of 0(m 4 ap ) or 
0(m 2 ap ,s). However, these additional terms have no effect on the final result, since only the 
unmodified s 2 behavior is extracted by the residue theorem. Because the tree-level amplitude 
Eq. (4.16) is analytic at the origin, we find \A"{s = 0) = — c c/q/| ao | 4 - The right-hand side 
of Eq. (4.18) involves an integral of a negative-definite quantity, as unitarity of the S'-matrix 
requires that Im[^4(s + ie)] < 0. 25 Hence, c must be positive. 26 

We conclude that the scalar self-interaction provides a negative contribution to the ef- 
fective potential at 0(g 4 T s /\ao\ 4 ). There are, of course, additional interactions from C^g 4 
that should also be considered. The local operators whose thermal expectation values lead to 
0(<7o^ 8 /l a o| 4 ) corrections to the effective potential are all dimension eight and involve four 
fields. They give rise to the set of two-loop diagrams shown in Figure 6. The scalar-scalar 
contribution I ss in Fig. 6 is the just-discussed I3. The other diagrams all have the same 
figure-eight topology and involve some pairing of scalars, fermions, and vectors running in 
the two loops. They were calculated in Ref. [3], and were found to all have the same relative 
sign. However, the sign of these contributions asserted in Ref. [3] is opposite to our conclusion 
for is S , stemming from the fact that the single overall coefficient c is negative in Ref. [3]. 2 
We now review the calculations of these diagrams and adjust the conclusions in lieu of the 
fact that c is actually positive. 

The F 4 terms of the low energy effective Lagrange density are readily obtained from the 

25 One may phrase the fact that the imaginary part of the forward amplitude must be negative in terms of the 
optical theorem, as done in Ref. [4] . Ultimately, the negative sign reflects S'-matrix unitarity, since 5 = 1 — iM 
and = 1 imply that Im(M) = -§|-M| 2 < 0. See, for example, problem 17 in Ch. 3 of Ref. [22]. 

26 As discussed in Ref. [4], this constraint on the sign of c is a special case of a more general scenario. The 
forward amplitude A(s), away from the real axis and probing mg ap -C s -C M^r, has a Taylor expansion 
around any point so in this region that begins as (s — so) 2 with a coefficient that is negative, up to corrections 
which scale as O(|so| 2 /M^, m 2 ap /M^). In our case, the low energy theory is weakly-coupled and this permits 
an analytic expansion for A(s) at the origin. 

27 The negativity of c in Ref. [3] may be traced back to one of the references cited in that paper. The value 
of c may be obtained from Eq. (3.11) of Ref. [18]; expressing their equation in the form of our Eq. (2.6) implies 
that c = — l/(87r 2 ), which has the wrong sign. 
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last term in expression (2.9) involving four copies of the spinor-valued field strength. The 
purely gauge interactions that are suppressed by no more than |ao| -4 are 



1 



r n=4 _ c % 

lb do 



{F^F^f - {F^FfJ) 2 . (4.19) 



The corresponding diagram in Figure 6 is I gg . We find that (for an arbitrary Lorentz gauge- 
fixing parameter), 

l | «o I 

where fJ++ = |^ + fl + = |«/+ and d = A. This is the same double sum-integral that we 
found in our evaluation of the scalar self-interactions. Note that I gg = I ss . 

The four-fermion interaction with all ^'s is easily derived from the same superfield in- 
tegral used to obtain the purely scalar interactions, namely expression (4.10). The purely tp 
interactions that are suppressed by |ao| are 

CfT 4 ^ -fnr [(^d^)(d v ^) - (^)( W)] • ( 4 - 21 ) 

The four-fermion interaction with all A's must have exactly the same form by SU (2)r sym- 
metry. In Figure 6, the corresponding diagrams are 1^ and We find 

4 

Ui> = h\ = 2^pV , (4.22) 
\ao\ 

where ^ = and J_ = = -7^/240. 

It is straightforward, but tedious, to find the interactions that mix scalars and vectors, or 
mix different types of Weyl fermion. To save ourselves some trouble we rely on the component 
field expression of given in Eq. (3.6) of Ref. [3]. Based on our discussion in footnote 

27, we shall factor out an overall — l/(87r 2 ) and replace it by the constant c. According to 
Ref. [3], the scalar-gauge and if)-\ interactions that are suppressed by |ao| -4 and contribute 
to finite temperature effects are given by 28 

Cff 4 => [{d^){d v ~a)F^ p F vp - (A^)(^<UA) - (Xd^)(^d u X)) . (4.23) 

|&o I 

In Figure 6, the corresponding diagrams are I sg and We find 

4 4 

I sg = 2(d-2)^(3Vn ++ , I lpx = 4-^(3Vn__. (4.24) 



28 Our er-matrix conventions differ from those of Ref. [3]. In particular, amours) = —oe (theirs), with similar 
spatial matrices. The two conventions are related by a spatial parity transformation which has the effect of 
conjugating cr-matrices by cr° (or <r° as appropriate). We have done a parity transformation in order to write 
the mixed ip-X interactions in our convention. 
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(Reassuringly, the gauge dependence cancels completely in I sg .) 

It remains to compute the diagrams for the scalar-fermion and gauge-fermion interactions. 
In Figure 6, these contributions are I s ^, I s \, I g ^, and I g \. Since each vertex involves one type 
of Weyl fermion paired with its Hermitian conjugate, SU{2)r symmetry requires equivalent 
interactions for ip and A. It follows that = I s \ and 1^ = I g \. The latter relation 
means that it is impossible for there to be any nontrivial gauge dependence in the diagrams 
involving gauge fields since there are no other diagrams left to cancel it. Since the complex 
scalar and gauge fields belong to the same supersymmetry multiplet, they have equal numbers 
of propagating degrees of freedom. Hence, all four diagrams must be equal. To determine 
their common value consider computing the index tr ((— l) F e~ /3H ) at weak coupling. In 
perturbation theory, the 0(<7o) contribution to the index comes from the class of diagrams 
shown in Figure 6, but with periodic temporal boundary conditions for all fields. This means 
that all frequency sums are taken over even integer multiples of ttT. Effectively, this changes 

all instances of fi__ and fi_] to The index, which must be an integer, cannot change 

as the coupling go is varied. Therefore, the 0(g 4 ) part of the index must be identically zero, 



— [Igs + Igg + 1^ + I\\ + / S g + + hip + IbX + Iglj) + 

16^4^0 ++ + 4J # | pbc/ 



p.b.c. 



(4.25) 



I ao 

Consequently, 

^| p . b . c . = -4|^^0 ++ . (4.26) 

This is a contribution to the index, but what we really want is the contribution to the partition 
function tr (e~^ H ). The functional representation of the trace requires antiperiodic temporal 
boundary conditions for fermions. Since each diagram involves a single fermion loop, we only 
need to turn one of the frequency sums in Eq. (4.26) into a sum over odd integer multiples 
of 7rT. Thus, the thermal contributions are given by 

4 

hi, = Isx = hi, = I s x = -4 ^% pV 0+_ , (4.27) 

l°o| 

where «+_ = = }j + J_. 

Having deduced the values of the diagrams in Figure 6, and knowing that c is positive, 
we can now understand how N = 2 gauge theory equilibrates at low temperature. The free 
energy density, viewed as a functional of ao, is the effective scalar potential after having 
integrated out all thermal fluctuations. Adding the blackbody and (undetermined) higher- 
order contributions yields 



(F(a )/V") ncutra i 



" 2 ^< T \ C (^ T " 
12 24 |o | 4 V |a | 4 



(4.28) 



The 'neutral' subscript is just a reminder that this is the contribution from the neutral degrees 
of freedom described by the low-energy effective Abelian theory; the heavy charged degrees 
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of freedom add the Boltzmann suppressed contribution (4.1). A key feature of the result 
(4.28) is that the free energy density decreases (becomes more negative) as one moves toward 
smaller values of |ao|- 

The subleading — T 8 /|ao| 4 behavior of the free energy density has a three-loop origin 
in the microscopic description of the theory. Analogous behavior is also observed in IIB 
supergravity calculations for the semiclassical region of SU{2) N = 2 theory [23] and in the 
SU(N C ) AA = 4 theory with N c — > oo and strong 't Hooft coupling [24]. 29 It is also worth 
noting that the F A interaction perturbs the free Hamiltonian by 



which is negative semi-definite since c is positive. Therefore, the F terms lower the classical 
energy. [This does not mean that the spectrum is unbounded below — the signs of higher 
powers of F 2 become important for large field strengths.] Reassuringly, similar behavior is 
found in other effective theories with Abelian gauge fields (e.g., the Born-Infeld action for a 
U(l) gauge field localized to a D-brane, or the Euler-Heisenberg action for QED) [4]. 

It is also reassuring to note that the sum of the free energy density contributions at low 
temperature from charged and neutral fields is consistent with the high temperature result. 
More precisely, the expressions for F/V far out on the free energy surface at T <C Myy, and at 
T S> M\y, match to leading order in the coupling. This follows from applying the asymptotic 

formula h(p) 4(2p) 3/2 e _7rp (derived in Ref. [9]) to Eq. (A.17), and comparing that to 

the sum of Eqs. (4.1) and (4.28). 

Lastly, let us make explicit why instanton contributions to the prepotential may be 
ignored compared to the non-holomorphic function fC. Far out on moduli space, one-instanton 
corrections to J~(a) take the form fx a 2 (A/a) 4 with /i / [1]. The off-shell form of C^g 2 g iven 
in expression (2.14) involves only the real and imaginary parts of J~"(a) and its derivatives. 
After inserting Eq. (4.2), suitably redefining the coupling constant to include one-instanton 
effects, and expanding in powers of a/an, the terms containing f\ and at least one field a 
begin at 0(A 4 /a[j). This leads to a subleading power correction relative to the leading result 



4.2 Near the massless monopole singularity 

We now switch attention to the strongly-coupled region of moduli space. Vacua near u = uq 
have spectra that include two types of BPS states: electrically charged W bosons with mass 

29 In Ref. [24], the supergravity interaction potential between a stack of N c coincident non-extremal D3-branes 
and a single "probe" D3-brane was computed. This potential was interpreted as arising from the Wilsonian 
effective action for the massless modes obtained by integrating out the massive modes of M = 4 theory on its 
Coulomb branch. The potential was found to be attractive. This implies that the leading dependence of the 
free energy on the scalar expectation value also comes with a minus sign. Indeed, Eq. (3.6) of Ref. [24] shows 
this explicitly. It was also argued that the weak-coupling expansion for the free energy density contains a 
nontrivial T 8 /|ao| 4 term. 




(4.29) 



(4.28). 
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M\y/A ~ O(l) and magnetically charged monopoles with mass M m /A ~ 0((u/uq) — 1). We 
assume the temperature is far below the strong scale, T <C A, but this leaves the freedom 
to consider two distinct regimes: (i) T 3> M m (hot monopoles), or («) T <C M m (cold 
monopoles). 

In terms of the monopole dynamics, case (i) is a high temperature regime, so it is natural 
to construct a three-dimensional effective theory as in Appendix A. Recall that the low energy 
theory near u = uq is an Abelian gauge theory of Ad = {AdiWdo) with hypermultiplet 
matter 7i = (Q, Q') of mass M m = \[2 |od|. The monopole couples locally to the dual photon 
so this is simply an N = 2 generalization of QED in four dimensions. The effective theory is 
infrared free with a coupling go = go(M m ) <1. In M = 1 superspace, 



-g 2 D £ Q ED = (J d 2 0\W%W Da + n.c\ + J d 2 9d 2 6A ] D A D 

+ j d 2 9d 2 9 (Q^e 2VD Q + Q'^e- 2VD Q') + (-iV2 j d 2 6 Q'A D Q + H.c.^ . 



(4.30) 



The chiral multiplets Q and Q' are oppositely charged under the magnetic U(l) gauge group, 
and under an ordinary f7(l)f flavor symmetry. The superpotential is uniquely fixed by M = 2 
supersymmetry. Under U(1)r transformations, Ad has charge 2, Wo a has charge 1, and 
Q and Q' are neutral. It follows that both Weyl fermions from the hypermultiplet have R- 
charge —1, so conservation of the C/(l)i?-current is anomalous at the one-loop level. One may 
determine the residual symmetry from the fact that td ~ -Mn(ai)) must be 2-7r-periodic 
under shifts of the effective theta angle. The global symmetry is SU(2)r x x ?7(l)f- In 

components, 



9d £qed = jFj^F^ + iXoa^df.XD + iipD^d^tpD + |<9 M a D | 2 



+ 



2(q\Dlpq + q'*\Dlpq< - qipDlpq' ~ q'^D^q ~ dDlpqlpq') + H.C. 



(4.31) 



+ 2\a D \ 2 {\q\ 2 + \q'\ 2 ) + \(\q\ 2 + \q'\ 2 )\ 



where F^ = d^A D — d u A D and = iziAjj^. The mass term M m for the hypermultiplet 
components appears when ao attains a translationally invariant expectation value. 

This N = 2 QED theory is valid below the momentum scale A. The next most relevant 
energy scale is the temperature T. Integrating out thermal fluctuations produces a three- 
dimensional effective theory which we denote as "QED3." By construction it will reproduce 
gauge invariant correlators for distances large compared to T" 1 . It is given by 

Z= [ VA^VA^VaoVqVq 1 expf ]- j d 3 x£ QE D 3 l , (4.32) 

J 1 9D,3 JV J 
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with 

£qed 3 = / + \{F%) 2 + U^ A °d) 2 + K(A) 2 + \%a D \ 2 + m 2 \a D \ 2 

+ \D+q\ 2 + \D7qf + (m 2 + (A° D ) 2 + 2\a D \ 2 ) {\q\ 2 + \q'\ 2 ) + \{\q\ 2 + \q'\ 2 ) 2 (4.33) 

>aD,q,q )■ 

The construction is similar in spirit to that for ESYM discussed in Appendix A. The fields 
are all mass dimension one bosonic zero-frequency modes that have been rescaled so that 
their kinetic terms have canonical normalization. To leading order in the dual coupling go, 
9d 3 = 9r>T- The covariant derivative acting on charged scalar components (originally from 
the hypermultiplet) is Df = di ± iAd% and the field strength is = d l A 3 L 



D 



d^A D . The 



effective theory QED3 has U{\) gauge invariance, plus translation and rotation symmetry. 
The global symmetries are realized as follows: ( „/* ) transforms as a doublet of SU(2)r, 
an — > —an under the (2^)^ generator, and q — > e luJ q and q' — > e~ luJ q' under U{1){ for arbitrary 
real oj. Gauge invariance allows mass terms for the various three-dimensional scalar fields. 
i?-symmetry requires that the hypermultiplet scalars appear in the invariant combination 
\q\ 2 + W\ 2 i an d that there are no operators cubic in ap. All other local, gauge invariant 
operators of mass dimension 4 or higher are lumped into (^thermal- 

The electrostatic mass raE, hypermultiplet mass m^, and dual scalar mass m s are fixed 
via matching calculations. One finds 



ml 



ml 



2mi 



g 2 D T 2 + 0(g%T 2 ). 



(4.34) 



Terms in 5f7 t hermai m ay be calculated using background field methods. Note that the tree-level 
scalar potential in the QED Lagrange density, given by the last line of Eq. (4.31), vanishes 
when q = q' = 0, regardless of the value of ao- So in the four-dimensional action one may 
expand around a saddle point = aoo (constant) and all other fields zero. Integrating out 
Gaussian fluctuations around this background leads to the following effective potential from 
non-static modes 30 



{F/g 2 D;i)U thermal (<2Z50 
T 2 



7T 7T 

— T H 

6 4 



all other fields zero 



M 2 
m +ln2 



n 2 T 2 

where the effective monopole mass 



M 2 

m 

7T 2 T 2 



n=3 



M 2 

m 

7T 2 T 2 



T 4 + 0(^T 4 ), (4.35) 



Mi 



2|a_oo| 



(4.36) 



An expression for the coefficients c n is given in Appendix A. In [/thermal] the constant term 
represents the blackbody radiation from an Abelian vector and hypermultiplet, and the co- 
efficient of the term quadratic in a^o agrees with m 2 ; everything else constitutes <5C/ t hermai- 



The mean-field-dependent quadratic forms in the shifted Lagrange density involve only the hypermultiplet 
component fields, so one does not need to fix a gauge at leading order. 
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Consider the momentum hierarchy T 3> M m S> g^T. 31 To integrate out massive 
monopole fields in QED3, expand the dual scalar field around its expectation value as od = 
{old) + cr with (an) = cldo- Then 



£qed 3 = a D (-V 2 + ml) a* D + \A% (-V 2 + ml) A% 

'-V 2 + m 2 + M 2 




+ [q, q' 




(4.37) 



where the ellipsis indicates terms cubic and higher order in fluctuations. The static hyper- 
multiplet contribution to the effective potential is 



(T/g 2 Dt3 )U static = 4/(M 2 )T[1 + 0(m 2 /M 2 )] , 



(4.38) 



where the function / is given in Eq. (A. 14). The overall factor of 4 accounts for the four real 
degrees of freedom in q and q' . 

The new lowest energy effective theory, valid for distances large compared to M~ l , is a 
three-dimensional U(l) gauge theory with coupling g DS which also includes a neutral real 
scalar A D with mass the and a neutral complex scalar a£> with mass m s . The free energy 



1 D 

density is obtained from the sum of ^/thermal and C/ S t a ticj 
F(a D0 )/V = T 4 ' 



IT 2 7T 2 



Ml 

m +ln2 



n 2 T 2 



M 2 

m 

7T 2 T 2 



n=3 



M 2 

m 

TT 2 T 2 



+ 0(g 2 D/ 



+ MlT 



1 



3tt 

+ 0(( to T) 3 T). 



+ 0(g 2 D T 2 /Ml) 



(4.39) 



Each line in Eq. (4.39) displays a contribution from one of the three momentum scales: T, 
M m , and gjyT (in that order). Defining a dimensionless mass ratio, 



Mm 
TFT 



we have 



^(a^o)/^ 



~ + ~h{p) + 0{gl) 



(4.40) 



(4.41) 



where h(p) is a function that increases monotonically for all p and is given explicitly by 
Eq. (A. 18). In the interval 1 > p > gz>, where Eq. (4.41) is valid, the free energy density is 
minimized as the monopole mass approaches zero. 32 

Near the monopole singularity, ajjo is mapped back to the gauge invariant coordinate 
u via the linear relation aoo ~ cq(u — uq), where cq = i/(2A) may be determined from 



31 The other regime, M m <JC goT, is unremarkable since, for sufficiently small ao, the O(goT) screening 
mass provides a big curvature at the origin of field space. 

32 Ref. [3] claims that a nontrivial minimum of the free energy density exists at M m ~ 0(gr>T). Since 
M m oc \u — «o|, this would imply an entire circle of minima in the u- plane centered around the massless 
monopole point. A continuous set of degenerate equilibrium states suggests a spontaneously broken continuous 
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the elliptic curve solution [1]. Since the free energy density decreases as u approaches uq, 
the point uq at which monopoles become massless must be a local equilibrium state. The 
effective theory at u = uq is infrared free, so the free energy density at this particular point 
is simply the blackbody contributions from a massless vector multiplet and hypermultiplet, 
up to corrections suppressed by the strong coupling scale, 

F/V = -^T 4 (1 + 0(T/A)). (4.42) 

Using the discrete i?-symmetry, this formula for F/V must also hold at u = —uq, the point 
where dyons become massless. Hence, there are two degenerate local minima of the free 
energy surface. 

Finally, let us consider case (ii), where T <C M m <C A, so the monopoles are cold and 
heavy and must be integrated out before considering the effects of thermal fluctuations. The 
resulting effective theory is given to next-to-leading order in the derivative expansion by a 
Lagrange density 

£ A eff = £ n D =2 cS + Cl=t s + 0(n > 6). (4.43) 

This describes the interactions of a massless U(l) M = 2 vector multiplet Ad = {Ad,Wdo) 
for distances large compared to M" 1 . As in Sec. 4.1, one may expand around a translationally 
invariant background od(x) = ano, define the small coupling g 2 D0 = 87r 2 /(ln | A/cido\ 2 — 3), 
and compute the free energy density perturbatively in qdq and as an expansion in inverse 
powers of a do. Since the leading logs in J 7 (a) and Fd{<id) only differ in functional form 
by a multiplicative factor, it follows that operators from £ D = g S do not contribute to F/V 
(aside from the trivial blackbody terms) until possibly O(g D0 T s /lapo] 4 ) [3]- 33 The leading 
correction to F/V comes from jCJTcff- No new work is needed to find the correction because 
electric-magnetic duality in the form (2.22) implies that KLo and K, have identical formulas. 
Hence, we simply adapt the result from Eq. (4.28). The free energy density functional is 

T 4 . (4.44) 

This expression decreases as one moves toward the massless monopole point, and crosses over 
into the form (4.41) when the monopole mass drops below T. 

symmetry. What is it? There are only two possibilities: U(l) C SU(2)r or C/(l)f, but an is not charged under 
either group. The resolution to this puzzle lies in properly handling the contributions of bosonic zero frequency 
modes to the effective scalar potential. The assertion of Ref. [3] follows from balancing the positive 0(M m T 2 ) 
term in Eq. (4.39) against the negative 0(g D M ln T 3 ) term, and this seems to imply that a minimum occurs 
at Mm/T ~ O{go). However, this is outside the range of validity of the calculation and is a misuse of the 
renormalization group. Recall that Eq. (4.39) is valid in the regime 1 3> M m /T 3> go- For M m /T <C go, 
the Wilsonian approach requires integrating out fluctuations whose correlation lengths are set by the Debye 
screening length (gnT) -1 , not M" 1 . Therefore, it is impossible to generate effective interactions which are 
non-analytic in M^. In particular, there is no term linear in M m in the effective potential. 

33 Corrections to To of the form A 2 ^2^ =1 ^(au/A)™ arise from integrating out infinitely many massive 
BPS states [25]. This can lead to 0(T/A) corrections in the free energy. We assume a separation of scales 
T < Mm C< Aso that T/A is still smaller than (T/M m ) 4 . 



F(a D0 )/V 



7T 

12 



24 odo 



\a>D0\ 
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5. Mass-deformed 517(2) Af = 4 gauge theory 

A simple generalization of pure N = 2 gauge theory is the addition of a single massive elemen- 
tary hypermultiplet in the adjoint representation. This theory (often referred to as N = 2*) is 
controlled in the far UV by a fixed point (the conformal N = 4 gauge theory) , but the relevant 
mass deformation induces running in the coupling, so that in the deep IR the theory is again 
pure M = 2 gauge theory. The Lagrange density can be obtained by adding to Eq. (2.2) the 
Kahler term for the hypermultiplet and the superpotential 

W = -i~ir(V2$[Q,Q']+mQQ'). (5.1) 

By a field redefinition, m may be chosen real. One is free to specify the value of an exactly 
marginal coupling in the UV. Let qo = e 2mT °, where tq = 9q/(2tt) + iAit/g 2 . is any complex 
number in the upper half plane. A choice of qo defines a scale Ao ~ |<?q m| where the 
theory evolves into the strongly-coupled pure TV = 2 theory [2] . We consider the limit of weak 
coupling, | go | *C 1, so that a large hierarchy exists between m and Ao- 

Classically, moduli space is given by Q = Q' = and [<£, = 0, so once again, vacua 
may be described as points in the -u-plane. Far from the origin, at \u\ 3> Aq, the weak coupling 
permits a mean field analysis. After applying the Higgs mechanism (for <f> = a<7 3 /2), one may 
read off the spectrum from the hypermultiplet F-term contributions to the scalar potential 
and the Higgs kinetic term. In N = 2 language, the vector multiplet splits into a massless 
photon .A 3 plus charged W bosons with masses \/2|a|. The hypermultiplet splits into a 
neutral component 7-L 3 with mass m and charged components 'Hr sl with masses |m=p y/2a\. A 
novel feature of this spectrum is that one of the electrically charged hypermultiplets (call it an 
"electron") can become massless at a = im/s/2. Therefore, in addition to the singularities 
where either a magnetic monopole or dyon goes massless, there is an additional singularity 
where an electron becomes massless [2]. This third singularity is located at u ~ \m 2 . A simple 
QED-like effective theory can be constructed near this point valid on distances 3> m~ 1 . By 
matching the SU (2) gauge coupling onto the one-loop QED coupling at the scale M\y ~ m, 
then running the QED coupling down to the mass scale of the light electron, one obtains the 
prepotential for an effective Abelian theory [26], 

Ha) ~ |r a 2 + ± a 2 Ing) -!-(«- m^f m( ^^g ) . (5.2) 

The perturbative analysis of Sec. 4.1 may be repeated for the prepotential (5.2) by ex- 
panding around a point near the third singularity, ao = m/ v2 + Aao . A similar conclu- 
sion is reached: the prepotential cannot contribute to the free energy density until at least 
O(f T*/\Aa n, where l/f = l/<? 2 + ^ln|^|. 

Instead of studying in detail four derivative terms in the effective action determined 
by a non-holomorphic function /C, one may argue that the third singularity must be a local 
minimum of the effective potential as follows. Since the electron becomes massless at«~ \m 2 , 
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the effective coupling go vanishes at this point and the low energy theory is free. For low 
temperatures, T -C Ao, the free energy density at the third singularity is simply F/V = 
_7T^ j>4^j _|_ 0(T/ Ao,T/m)). 34: Consider turning off the temperature and choosing a vacuum 
close to the third singularity. The spectrum still includes a massless photon, but the electron 
will have some small non-zero mass m e . Now turn on a temperature T <C m e . It is difficult 
to thermally excite electrons so their contribution to F/V will be exponentially suppressed. 
Since the effective theory near this point is weakly- coupled, we may trust the blackbody 
approximation to the free energy, but now this comes from a single type of N = 2 multiplet 
rather than two. Hence, the third singularity must lie deeper in the free energy surface than 
nearby points. 

Finally, we must understand the behavior of the free energy density for |u| 3> m 2 . One 
cannot simply adapt the J\f = 2 result since J\f = 2* is UV conformal, rather than asymptot- 
ically free. The non-holomorphic function fC for J\f = 4 theory is believed to be exactly of 
the form given in Eq. (4.9) without suffering renormalization [27, 28, 29]. The coefficient is 
known to be c = l/ir 2 [30]. With positive c, thermal fluctuations again make the asymptotic 
region of the free energy surface locally unstable. 

6. Discussion 

The non-trivial moduli space of SU (2) N = 2 supersymmetric Yang-Mills theory leads to a rich 
variety of dynamics on multiple length scales. Analyzing the thermodynamics of the theory 
requires careful application of effective field theory techniques to disentangle contributions 
from different types of fluctuations. At high temperature, T S> A, we found a unique Z2- 
invariant equilibrium state with a free energy density given by Eq. (3.4). At low temperatures, 
the flat zero-temperature ground state energy surface deforms into a non-trivial free energy 
surface. Far from the origin of moduli space, where Myy 3> A, we found that an arbitrarily 
small temperature, T <C Mw, causes the free energy surface to rise asymptotically. This 
corrects previously reported results [3] on the thermodynamics of this theory, and implies 
that minima of the free energy surface must lie in the portion of moduli space where the 
gauge coupling is strong. By using the dual description of the theory near the massless 
monopole (or dyon) points, we were able to analyze the thermodynamics when M m <C T <C A 
and monopoles (or dyons) are hot, as well as when T <C M m <C A and monopoles (or dyons) 
are cold. We found that the free energy surface has degenerate local minima at the massless 
monopole and dyon points. Our results are summarized in Figure 1(b). 

As there are no points in moduli space with enhanced gauge symmetry, the simplest 
scenario consistent with the above observations is to assume that there are no other local 
minima of the free energy density, and that the free energy surface smoothly decreases toward 
the massless monopole and dyon points throughout the intermediate regions in which no weak 
coupling description is applicable. This gives a simple, consistent picture in which the discrete 

34 Recall that free M = 2 vector and hypermultiplets each contribute —J2T 4 to the free energy density. 



- 31 - 



i?-symmetry is spontaneously broken at low temperatures, with two co-existing equilibrium 
states. 

The restoration of i?-symmetry at high temperature, combined with its spontaneous 
breakdown at low temperature, implies that there must be a genuine thermodynamic phase 
transition. The transition temperature must be some pure number times A, which is the only 
intrinsic scale in the theory. The spontaneous symmetry breaking, and consequent change 
in the number of distinct equilibrium states, ultimately arises from the existence of multiple 
special points in moduli space where equal numbers of massless states appear in the low 
energy spectrum. It is instructive to contrast this with SU(N C ) N = A gauge theory. The 
moduli space of this theory is locally flat and corresponds to the orbifold M. 6( - Nc ~ 1 ^ / Sn c . Only 
the single vacuum state at the origin is a fixed point of the entire permutation group. At 
this point the theory is superconformal and has the maximal number of massless gluons in its 
low energy spectrum. At weak coupling, these gluons provide the largest possible order one 
contribution to the free energy density in the form of blackbody radiation. At any non-zero 
temperature, there is a unique equilibrium state. 

We also examined weakly-coupled SU (2) M = 2* theory. At low temperature, we found 
that the additional hypermultiplet leads to the appearance of a third local minima in the free 
energy surface. This suggests the possibility of three distinct thermal equilibrium states. They 
correspond to points in the u-plane where a hypermultiplet (either solitonic or elementary) 
becomes massless. All three local minima on the free energy surface have the same free energy 
density, F/V = — ^ T 4 , up to corrections suppressed by the ratio of temperature to the strong 
coupling scale Ao- As only two local minima are related by the discrete i?-symmetry, we are 
unable to determine, based on our effective theory analysis, whether the three local minima 
are exactly degenerate, and if not which are the true global minima. 

An obvious extension of this work would be to generalize the low temperature analysis 
to SU(N C ) gauge groups with any number of colors. In particular, it would be interesting 
to study thermodynamics in the weakly-coupled SU(N C ) N = 2* theory. If an i?-symmetry 
phase transition exists in this theory, one could parametrize the transition temperature as 
T c = mf(N c , Ao/m) for some dimensionless function /. An understanding of the large N c 
limit of this formula might shed light on recent results obtained from the supergravity dual for 
the strongly-coupled version of the theory at finite temperature [31]. In the strong-coupling 
limit, the scale Ao is the same as the mass deformation m. Therefore, the critical temperature 
of the transition may be parametrized as T c = mf(N c ) for some unknown function /. A 
phase transition can be detected, in the large N c limit, by finding a zero of the free energy 
density as a function of m/T. However, numerical calculations do not show any zero-crossing 
behavior in the interval < m/T < 12 [32]. This is somewhat unexpected, as one might have 
expected qualitatively different behavior in the regimes T <C m and T>m [33] . One possible 
resolution would be for the function f(N c ) to scale as some (positive) power of 1/N C , forcing 
m/T c to move off to infinity as N c — > oo. Our work was originally motivated by the desire 
to better understand the nature of phase transitions in non-conformal gauge theories with 
supergravity duals. One of the basic questions for the thermodynamics of such theories is 
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understanding how the large N c limit affects the number and location of thermal equilibrium 
states. 
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A. High temperature effective theory 

In this appendix we compute the thermal effective potential for the complex scalar field 
(j) along its flat directions in J\f = 2 gauge theory. A unique global minimum is shown to be 
located at vanishing (j). Our procedure closely follows the treatment for weakly-coupled J\f = 4 
gauge theory in Ref. [9]. 

The most important effect of thermal fluctuations will be to generate temperature- 
dependent effective masses for certain fields. These are the complex scalar field (j) and the 
"time" component Aq of the real Euclidean gauge field. 35 At high temperature the effective 
coupling g{T) is small, so these masses are calculable in perturbation theory, are positive, and 
O(gT). 36 Consequently, there is positive curvature in field configuration space at (j) = Aq = 0. 
[This potential for Aq is consistent with the expected spontaneous breaking of center sym- 
metry at high temperature, since a vanishing mean and small fluctuations for Aq imply a 
non-zero mean for the Polyakov loop, tr V(e l f A ° dt ).] The point (j) = is a local minimum 
of the effective scalar potential and therefore also of the free energy surface. The goal is to 
now sharpen this observation and determine whether this local minimum is actually a global 
minimum of the effective scalar potential. 

Consider a regime where the temperature is much greater than the mass of W bosons, 
T S> M\y- Finite temperature effects on the long distance properties of the theory are 
captured by constructing a sequence of effective field theories. For an asymptotically free 
non-Abelian gauge theory at high temperature, the momentum scales shown in Eq. (3.1) 
are important for understanding its static properties. The temperature T is the typical 
momentum of particles in the plasma, the Debye scale gT sets the (inverse) correlation length 
of color-electric screening, and the magnetic mass gap g 2 T sets the (inverse) correlation length 

35 More precisely, one may view Ao as the traceless part of the gauge invariant Polyakov loop around the 
thermal circle [8]. 

36 Masses computed in strict perturbation theory should be understood as matching parameters in the low 
energy effective theory. The physical Debye screening mass, for example, includes nonperturbative corrections 
of order g 2 T [8]. 
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of color-magnetic screening. The separation between these scales is parametrically large at 
weak coupling. One may construct a low energy effective theory known as "ESYM" that 
reproduces static gauge invariant correlators of SU(2) M = 2 gauge theory for distances much 
greater than the inverse temperature [9]. 

The starting point for the construction is to consider the theory defined in a periodic 
spatial box of volume V and inverse temperature f3. The partition function has the functional 
integral representation 



Z = I VAn V\ £>V ^4> exp 



~ [ dxl [ d 3 xC 
9 Jo Jv 



(A.l) 



| 2 + ^ f l 2 



where C is given by the right hand side of Eq. (2.4). The Euclidean time direction is compact 
so the Hilbert space trace which defines the partition function requires periodic (antiperiodic) 
boundary conditions for bosonic (fermionic) fields. Each four-dimensional bosonic (fermionic) 
field may be decomposed as discrete Fourier series with frequencies that are even (odd) integer 
multiples of irT and coefficient functions that depend only on the three spatial directions. 
Integrating out all of the non-static modes (i.e., modes with frequencies \v\ > irT) results in 

Z = I VAiVA V<f> exp [- \ [ ^x^esymI , (A.2) 
J L 9i Jv J 

with 

£esym = / + 2tr Ui# + |(AA)) 2 + \vfi\A\ + |A0| 2 + m\\4>? + |[4>, 

+ ^thermal (-^j! A)> <?) ■ 

(A.3) 

This effective theory is a three-dimensional theory with SU(2) gauge invariance, translation 
and rotation symmetry, and Zg /^-symmetry which is realized as <fi — > e ln ^ 2 (p. 37 The Lagrange 
density (A.3) is obtained from dimensional reduction; this automatically produces the familiar 
gauge invariant derivative terms and commutator potentials. The adjoint covariant derivative 
Di = di + i[Ai, ■]. Gauge invariance allows mass terms for the three-dimensional scalar fields 
Ao and 4> — these terms are non-zero as a result of integrating out the infinite tower of non- 
static modes. The discrete i?-symmetry does not permit terms cubic in <j>. All other local, 
gauge invariant operators of mass dimension 4 and higher that can be constructed out of F%j , 
Aq, and (j) are contained in SUthevmai- The only operator with dimension zero is the identity 
and we include it explicitly with coefficient /. 

The three-dimensional fields of ESYM have engineering dimension one and represent 
renormalized zero-frequency modes of the corresponding four-dimensional fields. 38 One may 
choose their normalization so that £esym contains canonically normalized kinetic terms. 



37 SU(2)n invariance holds trivially since the fermion doublet has been integrated out, leaving only SU(2)r 
singlets. 

38 For example, the zero-frequency mode of the field Ao(x) is given by (x) = dx% Aq(x). The 
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To avoid unnecessary clutter in the presentation, we have chosen not to introduce separate 
notation for the renormalized fields of the effective theory. This slight notational sloppiness 
is harmless because the computation and results for the effective scalar potential presented 
in this appendix are unambiguous to leading order in the coupling. 

The coefficients of the operators in the effective theory are determined by computing 
correlators of static observables as a function of their spatial separation, at short distance, 
in the original N = 2 gauge theory and in ESYM. Perturbative calculations may be used 
to match results within the window 39 (irT)~ l <C \x— y\ <C (gT)^ 1 . The resulting ESYM 
effective theory, with all coefficients determined, then provides a description of physics valid 
for distances large compared to T _1 . 

The inverse 3d gauge coupling, 1/g 2 = Z 2 (3/g 2 , is an overall factor multiplying the action 
of the 3ci effective theory. We will only need to know g| to leading order, in which case the 
renormalization factor Z 2 may be replaced by unity. Hence g" 2 = g 2 T +0(g 4 T). Following the 
techniques described in the appendices of Ref. [9] , one finds that the coefficient of the identity 
operator, which represents the contribution to the free energy density from the momentum 
scales of order T and up [according to F/V = {T/g 2 )f + •••], is given by 



/ = 3<7 3 2 T 3 



(A.4) 



The leading term is just blackbody radiation while the 0(g 2 ) term comes from two- loop 
bubble diagrams. By computing the gluon self-energy tensor at zero external momentum in 
a loop expansion, one finds the electrostatic mass parameter wie. At leading order it is equal 
to the Debye mass, 

m 2 E = 2g 2 T 2 + 0(g 4 T 2 ). (A.5) 
A similar process determines the scalar mass parameter, 

m 2 =g 2 T 2 + 0{g i T 2 ). (A.6) 

The results (A.4)-(A.6) were also obtained in Ref. [12] by studying the thermodynamics 
of six-dimensional M = 1 gauge theory. This theory has eight supercharges and reproduces 
four-dimensional M = 2 gauge theory upon dimensional reduction. 



renormalized electrostatic field is defined by A T ell (x) = Z A ^A °\x), where Za is a dimensionless function of 
the coupling and fi/T, with fi an arbitrary renormalization scale. The wavefunction renormalization factor 
Za may be computed in perturbation theory and equals 1 + 0(g 2 \n(fi/T)). To the order at which we will be 
working this finite renormalization factor, as well as those of the other fields in the 3d effective theory, may be 
ignored. 

39 This range of scales is determined from the following considerations. Naive perturbative calculations in 
the short distance theory do not account for screening of color-electric fields in the plasma which cause Ao to 
develop a finite correlation length of order (gT) _1 . Far below this length scale, Ao behaves like a massless field 
and strict perturbation theory can be used [7] . The long distance theory includes the effects of integrating out 
thermal fluctuations with frequencies > ttT. Therefore, the effective theory is only valid for energies far below 
this scale. 
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Let us pause to understand our earlier claims about thermal masses for scalar fields. 
Like the zero temperature case, there is a commutator potential of the form [(ft, (ft] 2 that 
requires the expectation value of the three-dimensional field (ft to be a diagonal matrix (up to 
gauge transformations) . A new feature in the effective theory is the presence of the quadratic 
potential m^,|0| 2 . This term levies an energy cost for non-zero eigenvalues of (ft. For sufficiently 
small (ft, this term will completely dominate the behavior of the scalar potential. These 
tree-level considerations indicate that (ft = (along with Aq = 0) is a local minimum of the 
free energy density function. However, a tree-level analysis is not truly complete without 
understanding the effects of higher dimension operators in ^thermal- In principle one may fix 
the coefficients of these operators by matching successively higher n-point functions, but it is 
more convenient to use a background field method. 

The idea of the background field method is to choose a general constant value of the 
scalar field (other than (ft = 0) and evaluate the thermal effective potential as a function of (ft. 
We choose to expand (ft around (ft^t = ad 3 /2, with the gauge field and fermions vanishing, 
thus incurring no energy cost from gradients or tree- level potentials. Integrating out non- 
static fluctuations around this background field yields all higher order terms in the effective 
potential when the scalar field lies along its flat directions, 40 



{T/gi) ^thermal! 



9 flat : 



~-aa 3 /2) 



7T r^A 7T 
■— T + — 



all other fields zero 



\ 2 oo 

) +£< 



ra=3 



TT 



2j>2 



T 4 + 0( ff 2 T 4 ), (A.7) 



where M\y denotes the effective W mass, 



= 2\a\ 2 



(Ai 



and 



■ir i 



4-»)<S=8"c(2„-3). 



(2n)!! 



(A.9) 



The coefficients c n are derived in Appendix C of Ref. [9]. The quantity (T/g 2 ) [/thermal 
represents the contribution to the free energy density F/V from the momentum scale T (and 
above), and Eq. (A.7) is the result of integrating out only Gaussian fluctuations around the 
background field </>fl a t- In ^thermal) the constant term represents the blackbody radiation for a 
triplet of massless vector multiplets, and the coefficient of the term quadratic in a agrees with 
m|. We define 5C/ t hermai to be everything in [/thermal aside from the constant and quadratic 
terms. 



40 In the shifted action, there is a bilinear term arising from the covariant derivative of <ft that mixes the 
fluctuating part of <f> and A^. Integrating by parts shows that this term is proportional to the divergence of 
the gauge field. We choose Landau gauge, d^A^ — 0. Diagrammatically, the bilinear term generates a scalar- 
vector vertex with an external <^fl at leg at zero momentum. The Feynman rule for the vertex is proportional to 
the momentum of the fluctuating gauge field, but a Landau gauge vector propagator projects onto the space 
transverse to the momentum [13]. Thus, with this gauge choice one can ignore the bilinear mixing term. 
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At this point, we have only assumed that T S> Mjy. We now refine the restriction to 
T S> Myy S> gT. If a nontrivial minimum of the effective potential exists, then it must lie in 
this region. To see this, it is helpful to define a dimensionless rescaled mass, 

(-0) 

At lowest order in the coupling, expressed in terms of p, the contribution of non-static modes 
to the effective potential takes the form ^-/(/?)T 4 , where f(p) = — ^ + p 2 + (In 2) p 4 + 0(p 6 ) is 
completely independent of g. Including the effects of static modes will be shown to supplement 
/ with an additional cubic term, bp 3 . The coefficient b is negative and 0(g ). Therefore, a 
nontrivial minimum of the effective potential can exist only for some O(g ) value of p. 

To complete the analysis, one must consider the contribution to the free energy density 
from fluctuations of the static ESYM fields with masses set by M\y. One may do this by 
expanding cj) around its expectation value and integrating out all fields that develop finite 
correlation lengths of order Mw 1 . Let <fi = ((f)) + a with (<f>) = aa 3 /2 and a a fluctuating 
adjoint representation complex scalar field. This field redefinition naturally generates mass 
terms that depend on ((f)), as well as new interactions involving three or more fluctuating fields. 
Since the lowest order contributions to the free energy density are functional determinants 
obtained by integrating Gaussian fluctuations, we restrict our discussion to the quadratic 
forms in £esym obtained after the background field shift. For static gauge fields, 

£esym ^ \Al^A] + + M^W*, (A.ll) 

where Ajj = —5ijX7 2 + (1 — a~ 1 )didj is the covariance operator for vector fields and W{ = 
(A\ — iAj)/\/2. A Lorentz gauge-fixing term has been included. We choose to compute in 
Landau gauge (i.e., we send a — > + ). Expression (A.ll) indicates that an Mjy-dependent 
correlation length is generated for the off-diagonal components of the static gauge field. For 
the electrostatic scalar, 

£esym D ^o(-V 2 + ml)Al + W (-V 2 + m| + Af&)W *, (A.12) 

where Wq = (Aq — iAq)/\^2. Expression (A.12) shows that off-diagonal components of the 
electrostatic scalar also obtain Miy-dependent masses. For the complex scalar, 

£esym d , 3( _ v * + ml) ?> + E» (" V2 + <+ H 2 _ y2 -J + |o|2 ) E, (A.13) 



where £ = ^yji* Yo 2 *)' Since the mass matrix in expression (A.13) has eigenvalues 
and m 2 + M^, there is only one (complex) component of £ that obtains an M^-dependent 
mass. Notice that only off-diagonal components of the scalar field receive a-dependent mass 
shifts. 

Integrating over the off-diagonal fields in expressions (A.11)-(A.13) with correlation 
lengths of order M^ 1 is straightforward to leading order, since the integrals are Gaussian. 
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The basic ingredient needed to evaluate the resulting functional determinants is the single 
loop integral 

/( m 2) = i f i n (p + m 2 ) = -— (m 2 ) 3 / 2 . (A.14) 



(2nY 



12vr 



(We have regulated the theory by dimensional continuation to d = 3— 2e dimensions. 41 ) The 
contribution to the effective potential from the static modes with masses of order My/ is 
simply 

(T/gl) U static = 2 • 4 1(M^) r[l + 0( 5 2 T 2 /M^)] . (A.15) 

The factor of 2 accounts for the complex nature of all the off-diagonal fields. The factor of 
4 counts the total number of propagating fields: 2 for the gauge field since there are two 
transverse directions to a given spatial momentum, and 2 for the two kinds of scalars. The 
relative 0(g 2 T 2 /Myy) corrections in Eq. (A.15) come from evaluating / with the complete 
mass, either m 2 + Myy or + M^, and then expanding the result in powers of the small 
ratio (m^/My/) 2 or {m-^/My/) 2 ■ The minus sign in Eq. (A.14) is physically significant, as 
it shows that fluctuations of static fields have a destabilizing effect on the potential. If the 
magnitude of such destabilizing terms are large enough, then a nontrivial minimum will be 
produced away from the origin. 

The final low energy effective theory, valid for distances large compared to M^} 1 , is a 
three-dimensional U(l) gauge theory with coupling g 2 , a neutral real scalar Aq with mass 
die, and a neutral complex scalar ^> 3 with mass uia. The effective scalar potential in this 
theory, for (p along its flat directions, is given by the sum of ^/thermal and C/ S t a ti c - It follows 
that the free energy density, viewed as a functional of a, is given by 



F(a)/V 



7T 2 IT 2 

■T + y 



Mi 



w 



n 2 T 2 



+ ln2 



M^T 



3vr +U \ M? 



w / J 



+ 0({gTfT). 



7T 2 T 2 



2 oo 
+ J>„ 

n=3 



7T 2 T 2 



+ 0(g 2 



(A.16) 



Each line in expression (A.16) displays the contribution from one of the three momen- 
tum scales: T, Myy, and gT (in that order). The contribution from the soft scale gT 
may be obtained by integrating out the neutral scalars. In terms of the normalized mass 
p = y/2 \ a\/(irT), the final result is 



F(a)/V 



rp4 



(A.17) 



41 More precisely, one should first subtract from the integrand a term of the form ln(fc 2 + fJ. 2 ) where fi is an 
IR regulator. Alternatively, one may define I(m 2 ) by differentiating the formal integral with respect to m 2 . 
The resulting integral is only linearly sensitive to the UV cutoff. Analytically continuing in dimension and 
then integrating back with respect to m 2 leads to the stated result. 
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where 

h{p) = - 1 + p 2 - 4 p 3 + (In 2) p 4 + £ c n p 2 ". (A.18) 

n>3 

The function /i increases monotonically for all p (see Ref. [9] for a discussion of its global 
behavior). Our result for F(a)/V, which may be trusted in the interval 1 3> p S> g, is 
minimized as p approaches zero. 

In summary, the free energy density is given by Eqs. (A. 17)— (A.18) for 1 >^)>g. Since 
this function is monotonic, the free energy density is minimized precisely where the scalar 
field eigenvalues vanish. Hence, a = is the unique equilibrium state in the high temperature 
regime, T 3> A. 
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